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POLARIMBTRIC  PASSIVE  REMOTE  SENSING  OF  OCEAN  SURFACE 


Under  the  sponsorship  of  the  ONR  Contract  Contract  N0G014-92--J-1616  we  have 
published  4  refereed  journal  and  conference  papers. 

Update  of  progress  in  “Polarimetric  Passive  Remote  Sensing  of  Ocean  Surfaces” 

Our  work  on  this  project  in  the  past  year  has  concentrated  primarily  on  the  phase 
one  goals  as  outlined  in  the  original  proposal.  To  date,  we  have  performed  an  experiment 
in  which  the  polarimetric  thermal  emission  from  a  sinusoidal  water  surface  was  measured 
and  compared  against  a  theoretical  model,  performed  a  study  of  the  polarimetric  thermal 
emission  from  surfaces  randomly  rough  in  one  direction,  and  begun  work  on  a  model 
for  predicting  the  polarimetric  thermal  emission  from  surfaces  randomly  rough  in  two 
directions.  Each  of  these  studies  is  described  in  more  detail  below. 

Polarimetric  Thermal  Emission  from  a  Sinusoidal  Water  Surface 

To  verify  further  our  earlier  findings  reported  in  [1,2]  for  ocean-like  surfaces,  an 
experiment  was  designed  and  carried  out  in  which  the  polarimetric  thermal  emission  from 
a  sinusoidal  water  surface  was  measured  [3,4].  This  sinusoidal  water  surface  was  created 
by  placing  a  thin  sheet  of  fiberglass  with  a  sinusoidal  profile  on  top  of  a  pool  of  water 
and  removing  the  air  trapped  underneath  the  fiberglass  layer.  The  resulting  surface  was 
actually  a  ‘two-layer’  periodic  surface  whose  thermal  emission  should  be  close  to  that  of  a 
true  sinusoidal  water  surface  if  the  effect  of  the  fiberglass  layer  is  neglected.  The  first  three 
Stokes  parameters  of  the  thermal  emission  were  measured  at  both  10  GHz  and  14  GHz 
using  a  linearly  polarized  radiometer  whose  po’arzation  basis  was  rotated  to  perform  the 
polarimetric  measurements.  Significant  values  of  the  brightness  temperature  corresponding 
to  the  third  Stokes  parameter  U  were  observed  at  various  polar  and  azimuthal  angles  (as 
high  as  40  K  for  certain  configurations).  A  theoretical  model  for  the  thermal  emission  from 
such  a  two-layer  periodic  surface  was  constructed,  and  its  predictions  agreed  well  with  the 
expermental  measurements.  This  theoretical  model  also  indicated  that  the  fiberglass  layer 
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did  slightly  affect  the  brightness  temperatures  in  horizontal  and  vertical  polarizations,  but 
had  a  much  smaller  effect  on  U . 

Polarimetric  Thermal  Emission  from  Surfaces  Randomly  Rough  in  One  Direction 

Our  experiment  indicated  that  appreciable  values  of  U  could  be  obtained  from  a 
periodic  sinusoidal  water  surface.  However,  an  actual  wind  generated  ocean  surface  has 
a  very  complicated  structure  and  can  only  be  described  statistically  as  a  random  process. 
Thus,  a  numerical  experiment  in  which  the  polarimetric  thermal  emission  from  randomly 
rough  surfaces  was  investigated  was  next  performed  [5]. 

Many  approximate  theories,  such  as  the  Kirchhoff  approximation  and  the  small 
perturbation  method,  exist  for  use  in  predicting  the  thermal  emission  from  randomly  rough 
surfaces.  However,  the  theoretical  predictions  of  [1 J  indicate  that  large  surface  slopes  are 
needed  to  obtain  appreciable  values  of  U .  Since  the  approximate  theories  mentioned  above 
are  known  to  fail  if  the  surface  slopes  become  too  large,  a  Monte  Carlo  approach  was  chosen 
for  the  study. 

In  this  method,  a  set  of  surface  profiles  in  two  dimensions  with  given  statistics  was 
generated  using  a  spectral  method.  The  profiles  obtained  were  extended  infinitely  in  the 
third  direction.  The  thermal  emission  from  each  surface  of  the  set  was  then  calculated  using 
an  exact  integral  equation  approach.  The  results  from  the  set  were  averaged  to  obtain  the 
Monte  Carlo  estimate  of  the  polarimetric  thermal  emission.  The  surfaces  generated  were 
intended  to  model  the  ocean  in  the  microwave  frequency  range.  For  this  purpose,  a  power 
law'  spectrum  was  chosen  for  the  surface  statistics.  The  surface  profiles  generated  were  also 
made  periodic  to  avoid  problems  with  edge  conditions  in  the  calculations.  Convergence  of 
the  results  with  both  the  period  of  the  surface  chosen  and  the  number  of  surface  profiles 
averaged  was  investigated  in  the  study.  The  results  of  the  study  indicate  that  U  also  exists 
for  randomly  rough  surfaces  and  indicates  the  direction  of  the  surface  periodicity.  U  was 
also  found  to  be  sensitive  to  the  rms  height  of  the  surface  and  power  law’  spectrum  slope, 
and  insensitive  to  the  polar  angle,  permittivity  of  the  surface  medium,  and  high  frequency 
content  of  the  surface  spectra  beyond  a  certain  cutoff. 


Poiarimetric  Thermal  Emission  from  Surfaces  Randomly  Rough  in  Two  Directions 

The  preceeding  Monte  Carlo  study  gave  some  indications  as  to  the  properties  that 
the  third  Stokes  brightness  parameter,  U,  could  be  expected  to  possess  for  a  randomly 
rough  surface.  However,  the  effect  of  a  two  directional  roughness  spectrum  on  the  U 
predictions  is  not  yet  known.  For  this  purpose,  we  are  currently  working  on  a  model  for 
calculating  the  poiarimetric  thermal  emission  for  a  surface  rough  in  two  directions.  A 
Monte  Carlo  study  using  this  model  is  also  proposed  to  ascertain  the  properties  of  the  U 
parameter  for  this  more  accurate  ocean  surface  model.  This  model  is  based  on  modification 
of  the  extended  boundary  condition  approach  that  was  used  for  the  previous  Monte  Carlo 
study,  and  the  surface  statistics  to  be  used  will  model  both  the  upwind  and  crosswind 
components  of  the  ocean  wave  spectrum. 

Summary 

The  success  of  these  studies  in  measuring  and  theoretically  predicting  the  properties 
of  the  third  brightness  parameter,  U ,  have  advanced  the  idea  of  using  U  as  an  indicator  of 
ocean  wind  direction  from  remotely  sensed  data.  Our  work  in  the  development  of  a  model 
and  algorithm  for  the  realization  of  such  a  goal  will  continue  as  planned. 
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Recent  theoretical  works  have  suggested  the  potential  of  passive  polarimetry  in  the 
remote  sensing  of  geophysical  media.  It  has  been  shown  the  third  Stokes  parameter  U  of 
the  thermal  emission  may  become  large  for  azimuthally  asymmetric  fields  of  observation.  In 
order  to  investigate  the  potential  applicability  of  passive  polarimetry  to  the  remote  sensing 
of  ocean  surface,  measurements  of  the  polarimetric  thermal  emission  from  a  sinusoidal 
water  surface  and  a  numerical  study  of  the  polarimetric  thermal  emission  from  randomly 
rough  ocean  surfaces  were  performed. 

Measurements  of  sinusoidal  water  surface  thermal  emission  were  performed  using  a 
sinusoidal  water  surface  which  was  created  by  placing  a  thin  sheet  of  fiberglass  with  a 
sinusoidal  profile  in  two  dimensions  extended  infinitely  in  the  third  dimension  onto  a  water 
surface.  The  theory  of  thermal  emission  from  a  “two-layer”  periodic  surface  is  derived  and 
the  exact  solution  is  performed  using  both  the  extended  boundary  condition  method  (EBCJ) 
and  the  method  of  moments  (MOM).  The  theoretical  predictions  are  found  to  be  in  good 
agreement  with  the  experimental  results  once  the  effects  of  the  radiometer  antenna  pattern 
are  included  and  the  contribution  of  background  noise  to  the  measurements  is  modeled. 
The  experimental  results  show  that  the  U  parameter  indicates  the  direction  of  periodicity 
of  the  water  surface  and  can  approach  values  of  up  to  30  K  for  the  surface  observed. 

Next,  a  numerical  study  of  polarimetric  thermal  emission  from  randomly  rough  surfaces 
was  performed.  A  Monte  Carlo  technique  utilizing  an  exact  method  for  calculating  thermal 
emission  was  chosen  for  the  study  to  avoid  any  of  the  limitations  of  the  commonly  used 
approximate  methods  in  rough  surface  scattering.  In  this  Monte  Carlo  technique,  a  set  of 
finite  rough  surface  profiles  in  two  dimensions  with  desired  statistics  was  generated  and 
extended  periodically.  The  polarimetric  thermal  emission  from  each  surface  of  the  set  was 
then  calculated  using  both  the  EBC  and  the  MOM  and  the  results  were  averaged.  The 
surface  statistics  chosen  were  intended  to  model  a  wind  perturbed  ocean  surface  in  the  X 
to  Ku  band  microwave  region. 

The  results  of  the  study  indicate  that  the  U  parameter  is  sensitive  to  the  azimuthal  angle 
between  the  surface  periodicity  and  the  looking  angle  and  to  the  rms  height  of  the  surface, 
and  that  the  U  parameter  is  fairly  insensitive  to  variations  in  polar  angle,  permittivity, 
surface  power  law  spectrum,  and  surface  spectrum  high  frequency  cutoff.  These  properties 
give  further  strength  to  the  idea  of  using  the  U  parameter  to  detect  wdnd  direction  over 
the  ocean. 
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Abstract  A  numerical  study  of  the  polarimetric  thermal  emission  from  ocean  surfaces 
randomly  rough  in  one  dimension  using  a  Monte  Carlo  technique  is  presented.  In  this 
study,  a  set  of  finite  length  surface  profiles  with  desired  statistics  was  generated  using 
a  spectral  method.  Each  surface  was  extended  periodically  to  create  an  infinite  rough 
surface,  and  the  thermal  emission  was  computed  using  the  extended  boundary  condition 
method  (EBC)  and  the  method  of  moments  (MOM).  The  results  from  the  set  of  surfaces 
were  then  averaged  to  obtain  the  Monte  Carlo  estimate  of  polarimetric  thermal  emission. 

The  surface  statistics  chosen  were  intended  to  model  a  wind  perturbed  ocean  surface  in 
the  X  to  Ku  band  microwave  region.  The  results  of  the  study  show  that  the  third  Stokes 
parameter,  £/#,  is  sensitive  to  the  azimuthal  angle  between  the  surface  periodicity  and  the 
looking  angle,  the  rms  height  of  the  surface,  and  the  surface  power  law  spectrum  slope, 
and  that  this  parameter  is  insensitive  to  variations  in  polar  angle,  permittivity,  and  surface  I 

spectrum  high  frequency  content  as  an  indicator  of  the  azimuthal  asymmetry  of  the  surface. 

1  Introduction 

Recent  theoretical  works  have  suggested  the  potential  of  passive  polarimetry  in  the  remote 
sensing  of  geophysical  media  [1—5].  It  has  been  shown  that  the  brightness  temperature  in 
the  third  Stokes  parameter,  Ub,  may  become  large  for  azimuthally  asymmetric  fields  of 
observation.  In  references  [4]  and  [5],  values  of  Ub  as  high  as  30  K  were  measured  from 
sinusoidal  water  surfaces  at  10  and  14.6  GHz.  In  order  to  investigate  the  potential  appli¬ 
cability  of  passive  polarimetry  to  the  remote  sensing  of  rough  ocean  surfaces,  a  numerical 
study  of  the  polarimetric  thermal  emission  from  ocean  surfaces  randomly  rough  in  one 
dimension  was  performed.  A  Monte  Carlo  technique  utilizing  an  exact  integral  equation 
method  for  calculating  thermal  emission  was  chosen  for  the  study. 

The  next  section  presents  a  brief  background  on  the  theory  of  polarimetry  and  describes 
the  passive  polarimetric  brightness  vector.  The  method  of  calculation  for  the  numerical 
experiment  is  discussed  in  Section  3  and  the  numerical  results  are  presented  in  Section  4. 
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2  Theory  of  Polarimetry 

i 

In  passive  polarimetry,  brightness  temperatures  corresponding  to  all  of  the  four  modified 
Stokes  parameters  are  measured.  The  brightness  temperature  Stokes  vector  is  defined  as 
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where  Eh  and  Ev  are  the  emitted  electric  fields  received  from  the  horizontal  and  vertical 
polarization  channels  of  the  radiometer,  rj  is  the  characteristic  impedance,  and  C  =  K/X2 
with  K  denoting  Boltzmaam’s  constant,  A  the  wavelength.  The  first  two  parameters  of  the 
brightness  temperature  Stokes  vector  correspond  to  the  powers  received  in  the  horizontal 
and  vertical  polarization  channels,  respectively.  The  third  and  fourth  parameters  corre¬ 
spond  to  the  complex  correlation  between  the  electric  fields  received  by  the  horizontal  and 
vertical  channels.  We  will  label  the  four  parameters  Tbk,  Tb„,  Ub,  and  Vs  respectively  in 
this  paper.  It  is  shown  in  [3J  that  the  third  and  fourth  Stokes  parameters  may  be  related 
to  the  brightness  temperature  in  a  45  degree  linearly  polarized  measurement  ( Tbi )  and  a 
right  hand  circularly  polarized  measurement  (TBr)  as  follows: 


Ub  =  2  TSi  —  Tbh  —  TBv  (2) 

VB  —  2  Tbt  —  TBh  —  TBv  (3) 


Thus,  to  calculate  all  four  parameters  of  the  brightness  temperature  Stokes  vector,  the 
brightness  temperatures  in  horizontal,  vertical,  45  degree  linear,  and  right  hand  circular 
polarizations  are  first  calculated,  and  the  above  equations  are  used  to  obtain  UB  and  VB. 

In  the  passive  remote  sensing  of  rough  surfaces,  the  parameter  that  is  actually  of 
Interest  is  the  emissivity,  which  relates  the  brightness  temperature  emitted  by  an  object 
to  its  actual  physical  temperature,  under  the  assumption  that  the  object  is  at  a  constant 
physical  temperature  and  that  the  emission  from  the  object  is  the  only  source  of  brightness: 


TBa  =  ea(ff,  <f>)Tphv.  (4) 

In  the  above  equation,  the  subscript  a  refers  to  the  polarization  of  the  brightness  tempera¬ 
ture,  6  to  the  polar  observation  angle,  and  <f>  to  the  azimuthal  observation  angle.  Through 
the  principles  of  energy  conservation  and  reciprocity,  Kirchhoff’s  law  relates  this  emissivity 
to  the  reflectivity  of  the  surface  [6]: 

ea(6,<f>)  =  l -ra(6,<f>)  (5) 

The  reflectivity  ra(0,  <£)  for  the  given  incident  polarization  a  is  defined  as  the  fraction 
of  the  power  incident  from  direction  ( 0 ,  4>)  that  is  rescattered  and  can  be  evaluated  by 
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integrating  the  bistatic  scattering  coefficient  <f>\  8',<f>')  over  all  scattering  angles  in  the 

upper  hemisphere  and  summing  the  results  of  both  orthogonal  scattering  polarizations. 

ra(8, 4>)  —  ^2  r/2de' sine'  (6) 

4tt  ^  Jo  Jo 

In  the  above  expression,  (8,  <f>)  and  (#',<£')  represent  the  incident  and  the  scattered  direc¬ 
tions,  respectively,  and  the  subscripts  a  and  b  represent  the  polarizations  of  the  incident 
and  the  scattered  waves,  respectively. 

Thus,  to  calculate  the  fully  polarimetric  emission  vector,  the  bistatic  scattering  co¬ 
efficient  for  each  of  four  polarizations  is  first  calculated  and  integrated  over  the  upper 
hemisphere  to  obtain  the  reflectivity  for  that  particular  polarization.  Multiplication  of  the 
corresponding  emissivity  by  the  physical  temperature  of  the  object  under  view  yields  the 
brightness  temperature  for  this  polarization.  The  fully  polarimetric  brightness  vector  is 
then  calculated  as  described  previously.  A  physical  temperature  of  300  K  was  assumed  for 
the  surfaces  in  this  experiment. 

3  Method  of  Calculation 

Scattering  from  randomly  rough  surfaces  has  been  studied  extensively  by  a  number  of 
researchers.  Numerous  approximate  techniques,  such  as  the  Kirchhoff  approximation  and 
the  small  perturbation  method,  exist  and  work  well  under  certain  restrictions  for  providing 
the  statistically  expected  values  of  scattered  fields  and  power  [6],  However,  it  was  found 
in  [2j  that  in  order  to  measure  large  Ub  values  from  a  sinusoidal  surface,  a  large  height 
to  period  ratio  was  required.  The  large  slopes  of  such  a  surface  fall  into  the  regions  of 
non-applicability  for  the  above  methods,  so  a  Monte  Carlo  method  was  chosen  for  the 
study. 

In  this  Monte  Carlo  method,  a  set  of  surface  profiles  with  given  statistics  was  generated, 
and  the  fully  polarimetric  brightness  vector  from  each  surfa.ee  was  calculated  using  an 
exact  integral  equation  approach.  The  results  for  the  set  were  then  averaged  to  obtain  the 
final  Monte  Carlo  estimate  of  the  average  brightness  vector  for  those  particular  surface 
parameters.  The  surface  generation  and  analysis  procedure  are  described  in  more  detail 
in  the  next  section. 

3.1  Random  Surface  Generation 

A  spectral  method,  described  in  [7-8],  w'as  used  to  generate  a  finite  sample  of  a  rough 
surface  with  desired  statistics.  In  this  method,  a  set  of  normally  distributed  (both  real  and 
imaginary  part)  Fourier  coefficients  is  generated  and  then  weighted  by  a  spectral  density 
function  in  the  frequency  domain.  The  resulting  Fourier  coefficients  are  then  transformed 
back  into  the  space  domain  using  an  inverse  FFT  algorithm  to  obtain  a  surface  profile 
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sampled  in  space.  The  surface  profiles  were  linearly  interpolated  between  the  specified 
points. 

Since  only  a  finite  surface  profile  could  be  generated,  an  infinite  surface  was  created  by 
extending  this  finite  surface  periodically.  The  Floquet  modes  obtained  from  the  periodic 
surface  are  thus  a  discrete  approximation  to  the  continuous  spectrum  of  a  truly  infinite 
random  rough  surface.  The  effect  of  this  periodic  extension  on  the  polarimetric  brightness 
temperature  was  investigated  by  comparing  the  results  from  a  set  of  surfaces  with  a  given 
period  with  the  results  from  a  set  of  surfaces  a  longer  period  and  the  same  frequency 
spectrum.  A  period  of  20  wavelengths  was  found  to  be  sufficient  to  provide  convergence 
of  the  Ug  results  to  within  0.3  K  for  longer  period  surfaces.  The  surfaces  were  continuous 
when  made  periodic  due  to  the  periodic  properties  of  the  FFT  series  from  which  they  were 
generated. 

The  surface  profiles  generated  were  rough  in  one  dimension  only  due  to  *he  complexity 
of  the  calculations  for  surfaces  rough  in  two  spatial  directions.  The  rough  profile  along  the 
x  direction  was  extended  infinitely  along  the  y  direction  in  the  calculations  (see  Figure  1). 
Th  ree  dimensional  angles  of  incidence  were  allowed,  however,  so  that  polarimetric  effects 
could  be  observed  from  the  resulting  ‘conical  diffraction’  problem.  While  this  model  is 
extremely  simple  when  compared  to  .n  actual  wind  perturbed  ocean  surface,  its  results 
should  give  some  indication  as  to  the  properties  of  Ug  for  general  rough  surface*. 

The  surface  statistics  chosen  for  the  calculations  were  intended  to  model  ocean  surfaces 
in  the  microwave  frequency  range.  For  this  purpose,  a  dielectric  constant  of  50  +  :30  was 
chosen  for  the  medium  and  a  power-law  spectral  density  function  for  the  surfaces  was  used. 
This  function,  4>[k)  is  described  by: 


m  = 


—-<!■<  I- 

^  -  c(7) 

otherwise 


where  k  is  the  spatial  wavenumber  of  the  surface  in  ra  is/m ,  s  is  the  slope  of  the  power  law 
spectrum,  and  kc  is  the  cutoff  wavenumber  of  the  surface  set.  The  effects  of  varying  the 
slope  of  the  spectrum  and  its  high  frequency  cutoff  were  investigated  in  the  experiment. 
The  surface  in  Figure  1  is  one  of  the  surfaces  generated  with  -3  slope  power  law  spectrum 
and  a  high  frequency  cutoff  of  kc  =  An/X. 

Each  surface  was  normalized  independently  of  the  others  to  a  specified  rms  surface 
he:gh- :  rms  heights  of  A/20,  A/25  and  A/30  were  chosen  to  investigate  the  effect  of  varying 
tne  surface  height  on  the  polarimetric  brightness  temperature.  The  20A  requirement  for 
the  surface  period  prevents  larger  rms  heights  from  being  studied  due  to  computational 
limitations.  Each  generated  surface  consisted  of  400  points  within  this  peri'>i,  so  that  the 
surface  profile  was  sampled  every  A/20  along  the  x  direction.  This  sampling  frequency 
is  well  above  all  of  the  Nyquist  frequencies  of  the  power  spectral  density  functions  used 
in  the  experiment.  The  surface  statistics  calculated  from  the  generated  surface  sets  were 
compared  to  their  desired  theoretical  values  and  found  to  be  in  good  agreement. 

The  results  from  ten  surface  profiles  were  averaged  for  all  of  the  points  in  the  exper- 
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iment.  The  convergence  of  the  results  with  the  number  of  surface  profiles  averaged  was 
investigated  by  comparing  the  results  from  independent  groups  of  ten  surfaces.  The  results 
for  Ub  were  found  to  be  within  0.2  K. 

3.2  Integral  Equation  Method 

Numerous  studies  have  been  made  of  scattering  from  a  periodic  surface.  The  Extended 
Boundary  Condition  (EBC)  method  [9]  is  one  of  the  most  efficient  methods,  but  has  been 
found  to  become  ill  conditioned  for  steep  surfaces.  The  Method  of  Moments  (MOM) 
for  a  periodic  surface,  described  in  [2],  is  another  solution  of  the  problem  that  does  not 
have  the  steep  surface  problems  of  the  EBC,  but  is  more  computationally  intensive.  The 
calculations  from  these  two  methods  were  compared  for  a  sample  surface  for  each  case 
in  the  experiment  and  found  to  give  similar  results  in  all  of  the  cases.  This  agreement 
indicates  that  the  surface  slope  limitations  of  the  EBC  were  not  exceeded,  and  this  more 
efneient  r.ethod  of  calculation  was  used  for  the  larger  surface  sets  of  the  experiment. 


4  Results  of  Calculations 

Figure  2  is  a  plot  of  the  average  UB  value  for  the  A/20  , A/25,  and  A/30  rms  height  surfaces 
as  a  lunction  of  azimuthal  angle  for  a  polar  angle  of  20  degrees.  A  power  law  slope  of  —3, 
cutoff  wavenumber  of  kc  =  4ir/A,  and  dielectric  constant  of  50  -f  i30  were  used.  Points 
were  calculated  every  15  deg  in  azimuthal  angle,  so  that  the  interpolating  lines  shown 
should  not  be  taken  to  be  exact.  It  is  seen  that  the  Ub  value  is  small  for  viewing  angles 
parallel  or  perpendicular  to  the  direction  of  periodicity  of  the  surface,  and  that  the  highest 
values  of  Ur.  are  obtained  at  azimuthal  angles  of  approximately  45  or  135  degrees.  The 
magnitude  of  UB  is  determined  by  the  rms  height  of  the  surface,  with  rougher  surfaces 
giving  higher  values.  The  small  values  of  UB  obtained  in  these  results  (2  K)  are  due  to 
t.'u:  convergence  requirements  for  the  Monte  Carlo  technique:  larger  Ub  values  could  be 
oDtaincc  by  analyzing  rougher  surfaces,  but  would  require  more  computer  time  than  was 
avauab.e  for  this  study.  Also,  data  collected  from  the  SSM/I  and  reported  by  Wentz  [iQj 
indicates  tnat  the  azimuthal  variation  in  TBh  and  TBv  for  real  wind  perturbed  ocean  surfaces 
at  both  3/  and  19  GHz  is  typically  less  than  3  K,  as  is  obtained  for  these  rms  heights.  The 
5  channel  is  seen  to  indicate  the  azimuthal  direction  of  the  surface  in  all  cases  where  it  is 
•’rvalue.  Note  that  negative  values  for  Ub  are  possible,  as  this  parameter  corresponds 
omy  to  a  correlation  between  real  observables.  Also,  the  symmetry  about  <p  —  90  deg  seen 
in  figure  2  and  the  following  figures  is  due  only  to  the  statistical  symmetry  of  the  surface 
about  this  angle;  individual  surfaces  are  not  symmetric  about  <*>  =  90  deg.  This  gives 
an  iii  'ication  that  the  average  of  the  ten  surface  profiles  is  converging  to  the  true  statistical 


figure  3  investigates  the  effects  of  polar  angle  on  the  UB.  Plotted  are  the  results 
it  polar  fngles  of  20  ,40,  and  50  degrees  for  the  A/20  rms  height,  —3  slope,  k,  —  4x/A, 
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50 -ft  30  dielectric  constant  surfaces.  The  results  are  surprising  similar,  demonstrating  that 
this  parameter  is  relatively  insensitive  to  polar  angle  as  an  indicator  of  surface  azimuthal 
direction.  These  larger  polar  angles  are  important  due  to  the  need  for  large  polar  angles 
in  satellite  borne  sensor  applications. 

Figure  4  illustrates  the  effect  of  permittivity  on  Ub ■  Results  for  polar  angle  20  degrees, 
rms  height  A/20  ,  —3  slope,  kc  =  47r/A,  and  dielectric  constants  of  40  4  *25,  50  4  *30,  and 
60  4  *35  are  plotted.  The  results  again  indicate  that  the  Ub  channel  is  insensitive  to  a 
change  in  the  dielectric  constant  around  50  +  *30. 

Figure  5  displays  the  variation  in  Ub  with  the  slope  of  the  power  law  spectrum,  s.  The 
results  for  A/20  rms  height  surfaces  with  power  law  slopes  of  -3, -2. 5,  and  -4.5  are  compared 
for  a  polar  angle  of  20  degree,  cutoff  wavenumber  kc  —  4tt/A,  and  dielectric  constant  of 
50  4-  i30.  The  results  show  that  the  Ub  parameter  is  larger  for  the  -2.5  slope  surface,  ft 

which  has  more  high  frequency  content  and  hence  larger  slopes,  and  that  U b  is  smaller  for 
the  smoother  -4.5  slope  surface.  However,  Ub  exists  for  all  three  cases  and  indicates  the 
direction  of  surface  periodicity  in  each  case. 

Finally,  Figure  6  investigates  the  effect  of  the  high  frequency  cutoff  wavenumber,  kc  on 
Ub-  Plotted  are  the  results  from  A/20  rms  height  —3  slope  surfaces  for  a  polar  angle  of  20 
degrees  and  dielectric  constant  50  4-  t'30.  Cutoff  wavenumbers  of  4r/A,  8a/ A,  and  16?r/A 
are  shown.  This  change  in  the  cutoff  wavenumber  is  seen  to  have  little  effect  on  the  Up 
results. 

5  Conclusions 

A  numerical  study  of  polarimet.  -,  thermal  emission  from  randomly  rough  ocean  surfaces 
using  a  Monte  Carlo  method  has  been  carried  out.  This  study  indicates  that  poianmetnc 
information  can  be  useful  in  the  remote  sensing  of  anisotropic  rough  ocean  surfaces.  One 
proposed  application  is  the  remote  sensing  of  wind  direction  over  the  ocean,  as  the  rough 
surface  created  on  the  ocean  has  an  anisotropic  structure.  Since  the  Ub  parameter  is 
shown  to  be  relatively  insensitive  to  the  polar  angle,  the  permittivity  of  the  medium,  and 
the  cutoff  wavenumber  of  the  surface  spectral  density  function  as  an  indicator  of  surface 
azimuthal  direction,  a  wind  direction  sensor  could  feasibly  operate  under  a  wide  variety  of 
these  conditions.  Further  research  into  this  area  and  into  the  use  of  a  more  realistic  ocean 
surface  model  will  continue. 
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Abstract 

A  numerical  study  of  the  polarimetric  thermal  emission  from  ocean  surfaces  randomly 
rough  in  one  dimension  using  a  Monte  Carlo  technique  is  presented.  In  this  study,  a  set  of 
finite  length  surface  profiles  with  desired  statistics  was  generated  using  a  spectral  method. 
Each  surface  wa s  extended  periodically  to  create  an  infinite  rough  surface,  and  the  thermal 
emission  was  computed  using  the  extended  boundary  condition  method  (EBC)  and  the 
method  of  moments  (MOM).  The  results  from  the  set  of  surfaces  were  then  averaged  to 
obtain  the  Monte  Carlo  estimate  of  polarimetric  thermal  emission.  The  surface  statistics 
chosen  were  intended  to  model  a  wind  perturbed  ocean  surface  in  the  X  to  Ku  band 
microwave  region.  The  results  of  the  study  show  that  the  third  Stokes  parameter,  Ub, 
is  sensitive  to  the  azimuthal  angle  between  the  surface  periodicity  and  the  looking  angle, 
the  rms  height  of  the  surface,  and  the  surface  power  law  spectrum  slope,  and  that  this 
parameter  is  insensitive  to  variations  in  polar  angle,  permittivity,  and  surface  spectrum 
high  frequency  content. 
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1  Introduction 


Recent  theoretical  works  have  suggested  the  potential  of  passive  polarimetry  in  the  remote 
sensing  of  geophysical  media  [1-5].  It  has  been  shown  that  the  brightness  temperature  in 
the  third  Stokes  parameter,  Ub,  may  become  large  for  azimuthally  asymmetric  fields  of 
observation.  In  references  [4]  and  [5],  values  of  Ub  as  high  as  30  K  were  measured  from 
sinusoidal  water  surfaces  at  10  and  14.6  GHz.  In  order  to  investigate  the  potential  appli¬ 
cability  of  passive  polarimetry  to  the  remote  sensing  of  rough  ocean  surfaces,  a  numerical 
study  of  the  polarimetric  thermal  emission  from  ocean  surfaces  randomly  rough  in  one 
dimension  was  performed.  A  Monte  Carlo  technique  utilizing  an  exact  integral  equation 
method  for  calculating  thermal  emission  was  chosen  for  the  study. 

The  next  section  presents  a  brief  background  on  the  theory  of  polarimetry  and  describes 
the  passive  polarimetric  brightness  vector.  The  method  of  calculation  for  the  numerical 
experiment  is  discussed  in  Section  3  and  its  convergence  investigated  in  Section  4.  Finally, 
the  numerical  results  are  presented  in  Section  5. 


2  Theory  of  Polarimetry 


In  passive  polarimetry,  brightness  temperatures  corresponding  to  all  of  the  four  modified 
Stokes  parameters  are  measured.  The  brightness  temperature  Stokes  vector  is  defined  as 
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where  Eh  and  Ev  are  the  emitted  electric  fields  received  from  the  horizontally  and  vertically 
polarized  channels  of  the  radiometer,  tj  is  the  characteristic  impedance,  and  C  =  K/X2 
with  K  denoting  Boltzmann’s  constant,  A  the  wavelength.  The  first  two  parameters  of  the 
brightness  temperature  Stokes  vector  correspond  to  the  received  powers  in  the  horizontally 
and  vertically  polarized  channels,  respectively.  The  third  and  fourth  parameters  correspond 
to  the  complex  correlation  between  the  electric  fields  received  by  the  horizontally  and 
vertically  polarized  channels.  We  will  label  the  four  parameters  Tjgh,  Tbv,  Ub ,  and  Vb 
respectively  in  this  paper.  It  is  shown  in  [3]  that  the  third  and  fourth  Stokes  parameters 
may  be  related  to  the  brightness  temperature  in  a  45  degree  linearly  polarized  measurement 
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(Tg/)  and  a  right  hand  circularly  polarized  measurement  (Tgr)  as  follows: 


(2) 

(3) 


Ub  =  2  Tbi  —  Tbh  —  Tbv 
Vb  =  2  Tgr  —  Tsh  —  Tbv 

Thus,  to  calculate  all  four  parameters  of  the  brightness  temperature  Stokes  vector,  the 
brightness  temperatures  in  horizontal,  vertical,  45  degree  linear,  and  right  hand  circular 
polarizations  are  first  calculated,  and  the  above  equations  are  used  to  obtain  Ub  and  Vb- 

In  the  passive  remote  sensing  of  rough  surfaces,  the  parameter  that  is  actually  of 
interest  is  the  emissivity,  which  relates  the  brightness  temperature  emitted  by  an  object 
to  its  actual  physical  temperature,  under  the  assumption  that  the  object  is  at  a  constant 
physical  temperature  and  that  the  emission  from  the  object  is  the  only  source  of  brightness: 

TBa  =  ea(0, 4>)Tphy,  (4) 

In  the  above  equation,  the  subscript  a  refers  to  the  polarization  of  the  brightness  tempera¬ 
ture,  8  to  the  polar  observation  angle,  and  <f>  to  the  azimuthal  observation  angle.  Through 
the  principles  of  energy  conservation  and  reciprocity,  Kirchhoff ’s  law  relates  this  emissivity 
to  the  reflectivity  of  the  surface  [6]: 

ea(<9,  <f>)  =  1  -  ra(8,  <j> )  (5) 

The  reflectivity  ra(6,4>)  for  the  given  incident  polarization  a  is  defined  as  the  fraction 
of  the  power  incident  from  direction  ( 0,<f> )  that  is  rescattered  and  can  be  evaluated  by 
integrating  the  bistatic  scattering  coefficient  7(»(0,  <f>\  O' ,  <f>')  over  all  scattering  angles  in  the 
upper  hemisphere  and  summing  the  results  of  both  orthogonal  scattering  polarizations. 

ra(#,<^)  =  —  "jL,J0  J  d<t>'  7^(0 ,  <p-  O' ,  <j>')  (6) 

In  the  above  expression,  ( 0,4> )  and  (0',<p')  represent  the  incident  and  the  scattered  direc¬ 
tions,  respectively,  and  the  subscripts  a  and  b  represent  the  polarizations  of  the  incident 
and  the  scattered  waves,  respectively. 

Thus,  to  calculate  the  fully  polarimetric  emission  vector,  the  bistatic  scattering  co¬ 
efficient  for  each  of  four  polarizations  is  first  calculated  and  integrated  over  the  upper 
hemisphere  to  obtain  the  reflectivity  for  that  particular  polarization.  Multiplication  of  the 
corresponding  emissivity  by  the  physical  temperature  of  the  object  under  dew  yields  the 
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brightness  temperature  for  this  polarization.  The  fully  polarimetric  brightness  vector  is 
then  calculated  as  described  previously.  A  physical  temperature  of  300  K  was  assumed  for 
the  surfaces  in  this  experiment. 

3  Method  of  Calculation 

Scattering  from  randomly  rough  surfaces  has  been  studied  extensively  by  a  number  of 
researchers.  Numerous  approximate  techniques,  such  as  the  Kirchhoff  approximation  and 
the  small  perturbation  method,  exist  and  work  well  under  certain  restrictions  for  providing 
the  statistically  expected  values  of  scattered  fields  and  power  [6].  However,  it  was  found 
in  [2]  that  in  order  to  measure  large  Ub  values  from  a  sinusoidal  surface,  a  large  height 
to  period  ratio  was  required.  The  large  slopes  of  such  a  surface  fall  into  the  regions  of 
non-applicability  for  the  above  methods,  so  a  Monte  Carlo  method  was  chosen  for  the 
stud}'. 

In  this  Monte  Carlo  method,  a  set  of  surface  profiles  with  given  statistics  was  generated, 
and  the  fully  polarimetric  brightness  vector  from  each  surface  was  calculated  using  an 
exact  integral  equation  approach.  The  results  for  the  set  were  then  averaged  to  obtain  the 
final  Monte  Carlo  estimate  of  the  average  brightness  vector  for  those  particular  surface 
parameters.  The  surface  generation  and  analysis  procedure  are  described  in  more  detail 
in  the  next  section. 

3.1  Random  surface  generation 

A  spectral  method,  described  in  [7-8],  was  used  to  generate  a  finite  sample  of  a  rough 
surface  with  desired  statistics.  In  this  method,  a  set  of  normally  distributed  (both  real  and 
imaginary  part)  Fourier  coefficients  is  generated  and  then  weighted  by  a  spectral  density 
function  in  the  frequency  domain.  This  spectral  density  function  corresponds  to  the  desired 
spectral  density  of  the  surface  set;  it  is  also  the  Fourier  transform  of  the  desired  correlation 
function  for  the  set.  The  resulting  Fourier  coefficients  are  then  transformed  back  into  the 
space  domain  using  an  inverse  FFT  algorithm  to  obtain  a  surface  profile  sampled  in  space. 
The  surface  profiles  were  linearly  interpolated  between  the  specified  points. 

Since  only  a  finite  surface  profile  could  be  generated,  an  infinite  surface  was  created  by 
extending  this  finite  surface  periodically.  The  Floquet  modes  obtained  from  the  periodic 
surface  are  thus  a  discrete  approximation  to  the  continuous  spectrum  of  a  truly  infinite 


random  rough  surface.  The  effect  of  this  periodic  extension  on  the  polarimetric  brightness 
temperature  was  investigated  by  comparing  the  results  from  a  set  of  surfaces  with  a  given 
period  with  the  results  from  a  set  of  surfaces  with  twice  that  period  and  the  same  frequency 
spectrum.  A  period  of  2GA  was  found  to  be  sufficient  to  provide  convergence  of  the  results 
for  longer  period  surfaces,  as  will  oe  discussed  in  Section  4.  The  surface  was  continuous 
when  made  periodic  due  to  the  periodic  properties  of  the  FFT  series  from  which  it  was 
generated. 

The  surface  profiles  generated  were  rough  in  one  direction  only  due  to  the  complexity 
of  the  calculations  for  surfaces  rough  in  two  spatial  directions.  The  rough  profile  along  the 
x  direction  was  extended  infinitely  along  the  y  direction  in  the  calculations  (see  Figure  1). 
Three  dimensional  angles  of  incidence  were  allowed,  however,  so  that  polarimetric  effects 
could  be  observed  from  the  resulting  ‘conical  diffraction1  problem.  While  this  model  is 
extremely  simple  when  compared  to  am  actual  wind  perturbed  ocean  surface,  its  results 
should  give  some  indication  as  to  the  properties  of  Ub  for  general  rough  surfaces. 

The  surface  statistics  chosen  for  the  calculations  were  intended  to  model  ocean  surfaces 
in  the  microwave  frequency  range.  For  this  purpose,  a  dielectric  constant  of  50  4- 130  was 
chosen  for  the  medium  and  a  power-law  spectral  density  function  for  the  surfaces  was  used. 
This  function,  <f>(k )  is  described  by: 

2*/S  X<t<kc 


otherwise 


where  k  is  the  spatial  wavenumber  of  the  surface  in  rads/m,  s  is  the  slope  of  the  power 
law  spectrum,  and  kc  is  the  cutoff  wavenumber  of  the  surface  set.  The  effects  of  varying 
the  slope  of  the  spectrum  and  its  cutoff  were  investigated  in  the  experiment.  The  surface 
in  Figure  1  is  one  of  the  generated  s  —  3  slope  power  law  spectra  surfaces  with  a  high 
frequency  cutoff  of  kc  =  4tt/X  and  period  of  5A. 

Each  surface  was  normalized  independently  of  the  others  to  a  specified  rms  surface 
height;  rms  heights  of  0.075A,  0.05A,  0.033A  were  chosen  to  investigate  the  effect  of  varying 
the  surface  height  on  the  polarimetric  brightness  temperature.  The  20A  requirement  for 
the  surface  period  prevents  larger  rms  heights  from  being  studied  due  to  computational 
limitations.  Each  generated  surface  consisted  of  400  points  within  this  period,  so  that  the 
surface  profile  was  sampled  every  A/20  along  the  x  direction.  This  sampling  frequency 
is  well  above  all  of  the  Nyquist  frequencies  of  the  power  spectral  density  functions  used 
in  the  experiment.  The  surface  statistics  calculated  from  the  generated  surface  sets  were 
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compared  to  tlieir  desired  theoretical  values  and  found  to  be  in  good  agreement. 

The  results  from  twenty  surface  profiles  were  averaged  for  all  of  the  0.075A  rms  height 
points  in  the  experiment.  The  convergence  of  the  results  with  the  number  of  surface 
profiles  averaged  was  investigated  by  comparing  the  results  from  independent  groups  of 
twenty  surfaces  for  one  case  of  the  experiment.  The  results  for  Ub  were  found  to  be  within 
0.5  K  for  the  0.075A  rms  height  surfaces  at  polar  angle  20  degrees,  as  is  discussed  in  Section 
4. 

3.2  Integral  Equation  Method 

Numerous  studies  have  been  made  of  scattering  from  a  periodic  surface.  The  Extended 
Boundary  Condition  (EBC)  method  [9]  is  one  of  the  most  efficient  methods,  but  has  been 
found  to  become  ill  conditioned  for  steep  surfaces.  The  Method  of  Moments  (MOM) 
for  a  periodic  surface,  described  in  [2],  is  another  solution  of  the  problem  that  does  not 
have  the  steep  surface  problems  of  the  EBC,  but  is  more  computationally  intensive.  The 
calculations  from  these  two  methods  were  compared  for  a  sample  surface  for  each  case 
in  the  experiment  and  found  to  give  similar  results  in  all  of  the  cases.  This  agreement 
indicates  that  the  surface  slope  limitations  of  the  EBC  were  not  exceeded,  and  this  more 
efficient  method  of  calculation  was  used  for  the  larger  surface  sets  of  the  experiment. 

4  Verification  of  Convergence 

Two  convergence  issues  are  important  for  this  Monte  Carlo  method:  the  convergence  of 
the  results  with  the  number  of  surfaces  averaged  and  the  convergence  of  the  results  with 
the  period  of  the  surfaces  analyzed.  Both  issues  are  investigated  in  this  section. 

The  convergence  of  the  results  with  the  number  of  surfaces  averaged  is  the  first  type  of 
convergence  that  must  be  investigated.  To  check  this  type,  independent  groups  of  surfaces 
were  analyzed  and  their  results  compared.  It  was  found  that  the  difference  obtained  in  the 
average  b b  values  from  ten  surface  profiles  were  within  0.3  K  for  all  of  the  cases  with  rms 
height  less  than  or  equal  to  0.05A.  For  the  rms  height  0.075A  case,  twenty  surface  profiles 
were  averaged,  and  a  maximum  Ub  difference  of  0.5  K  in  the  average  from  independent  sets 
of  twenty  surfaces  was  observed.  Since  this  uncertainty  is  much  less  than  the  azimuthal 
variation  of  Ub  obtained  in  the  study,  these  averages  were  considered  sufficient  for  the 
study. 
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The  second  type  of  convergence,  convergence  with  respect  to  the  period  of  the  surfaces 
analyzed,  was  also  investigated.  This  convergence  was  verified  by  comparing  the  average 
results  from  groups  of  ten  surfaces  which  had  different  periods  but  the  same  frequency 
spectra.  In  effect,  this  procedure  corresponds  to  sampling  the  frequency  spectrum  at  a 
greater  number  of  points  in  the  frequency  domain  as  the  period  is  allowed  to  increase. 
Longer  scale  variations  were  cutoff  as  the  period  was  increased.  Figure  2  displays  the 
change  in  polarimetric  brightness  temperature  as  the  period  was  increased  for  0.05A  rms 
height  surfaces  at  a  polar  angle  of  20  degrees  with  power  law  slope  s  =  3  and  cutoff 
wavenumber  47r/A.  Plotted  are  the  results  for  surfaces  with  periods  of  5A,  10A,  15A,  20A, 
and  25A. 

The  characteristic  “jumps”  of  periodic  surface  scattering  are  observed  in  all  of  these 
curves,  but  it  is  seen  that  the  effect  becomes  less  pronounced  as  the  period  is  increased. 
This  is  due  to  the  fact  that  a  greater  number  of  Floquet  modes  is  obtained  for  the  longer 
period  surfaces,  which  results  in  a  smaller  shift  in  total  power  as  one  mode  makes  the 
transition  from  propagating  to  non-propagating.  Figure  3  compares  the  results  for  the  20A 
and  25A  period  surfaces  alone.  It  is  seen  that  these  curves  are  fairly  smooth  in  azimuthal 
angle  and  that  they  agree  to  within  0.2  K  in  Ub  f°r  ah  of  the  angles  analyzed.  The 
results  from  period  20A  and  25A  surfaces  were  also  compared  for  a  set  of  twenty  rms  height 
0.075A  surfaces.  For  this  case,  a  maximum  difference  of  1.5  K  was  obtained  in  the  Ub 
results.  However,  the  Ub  results  are  also  larger  for  this  case,  so  that  the  uncertainty  is 
approximately  the  same  percentage  of  Ub  as  in  the  smaller  rms  height  cases.  A  period  of 
20A  was  chosen  for  all  of  the  surfaces  throughout  the  rest  of  the  experiment. 

5  Results  of  calculations 

Figure  4  is  a  plot  of  the  average  Ub  value  for  the  0.075A  ,  0.05A,  and  0.033A  rms  height 
surfaces  as  a  function  of  azimuthal  ingle  for  a  polar  angle  of  20  degrees.  A  power  law 
slope  of  s  —  3,  cutoff  wavenumber  of  kc  =  47t/A,  and  dielectric  constant  of  50  +  t30  were 
used  in  the  calculations.  Points  were  calculated  every  15  degrees  in  azimuthal  angle,  so 
that  the  interpolating  lines  shown  should  not  be  taken  to  be  exact.  It  is  seen  that  the  Ub 
value  is  small  for  viewing  angles  parallel  or  perpendicular  to  the  direction  of  periodicity 
of  the  surface,  and  that  the  highest  values  of  Ub  are  obtained  at  azimuthal  angles  of 
approximately  45  or  135  degrees.  The  magnitude  of  Ub  is  determined  by  the  rms  height  of 
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the  surface,  with  rougher  surfaces  giving  higher  values.  The  small  values  of  Ub  obtained  in 
these  results  (5  K)  are  due  to  the  convergence  requirements  for  the  Monte  Carlo  technique: 
larger  Ub  values  could  be  obtained  by  analyzing  rougher  surfaces,  but  would  require  more 
computer  time  than  was  available  for  this  study.  Also,  data  collected  from  the  SSM/I 
and  reported  by  Wentz  [10]  indicates  that  the  azimuthal  variation  in  Tbh  and  Tqv  for  real 
wind  perturbed  ocean  surfaces  at  both  37  and  19  GHz  is  typically  less  than  3  K,  as  is 
obtained  for  these  rms  heights.  The  Ub  channel  is  seen  to  indicate  the  azimuthal  direction 
of  the  surface  in  all  cases  where  it  is  observable.  Note  that  negative  values  for  Ub  arc 
possible,  as  this  parameter  corresponds  only  to  a  correlation  between  real  observables. 
Also,  the  symmetry  about  =  90  deg  seen  in  Figure  2  and  the  following  figures  is  due 
only  to  the  statistical  symmetry  of  the  surface  set  about  this  angle;  individual  surfaces  are 
not  symmetric  about  <p  —  90  deg.  This  gives  an  indication  that  the  average  of  the  surface 
profiles  is  converging  to  the  true  statistical  average. 

Figure  5  investigates  the  effects  of  polar  angle  on  the  Ub-  Plotted  are  the  results 
at  polar  angles  of  20,  40,  and  50  degrees  for  the  0.075A  rms  height,  s  —  3,  kc  =  4rr/A, 
50-H30  dielectric  constant  surfaces.  The  results  are  surprising  similar,  demonstrating  that 
this  parameter  is  relatively  insensitive  to  polar  angle  as  an  indicator  of  surface  azimuthal 
direction.  These  larger  polar  angles  are  important  due  to  the  need  for  large  polar  angles 
in  satellite  borne  sensor  applications. 

Figure  6  illustrates  the  effect  of  permittivity  on  Ub-  Results  for  polar  angle  20  degrees, 
rms  height  0.075A  ,  s  —  3,  kc  =  47r/A,  and  dielectric  constants  of  40  -f  *25,  50  +  z'30,  and 
60  +  *35  are  plotted.  The  results  again  indicate  that  the  Ub  channel  is  insensitive  to  a 
change  in  the  dielectric  constant  around  50  +  *30. 

Figure  7  displays  the  variation  in  Ub  with  the  slope  of  the  powTer  law  spectrum,  s.  The 
results  for  0.075A  rms  height  surfaces  with  power  law  slopes  of  s  =  3,  s  =  2.5,  and  s  =  3.5 
are  compared  for  a  polar  angle  of  20  degree,  cutoff  wavenumber  kc  =  47r/A,  and  dielectric 
constant  of  50  +  i'30.  The  results  show  that  the  Ub  parameter  is  larger  for  the  ;  =  2.5 
surface,  which  has  more  high  frequency  content  and  hence  larger  slopes,  and  that  Ub  is 
smaller  for  the  smoother  s  =  3.5  slope  surface.  However,  Ub  exists  for  all  three  cases  and 
indicates  the  direction  of  surface  periodicity  in  each  case. 

Finally,  Figure  8  investigates  the  effect  of  the  high  frequency  cutoff  wavenumber,  kc  on 
UB-  Plotted  are  the  results  from  0.075A  rms  height  s  =  3  surfaces  for  a  polar  angle  of  20 
degrees  and  dielectric  constant  50  4-  *30.  Cutoff  wavenumbers  of  47r/A,  8tt/A,  and  16tt/A 


are  shown.  This  change  in  the  cutoff  wavenumber  is  seen  to  have  little  effect  on  the  Ujj 
results. 

6  Conclusions 

A  numerical  study  of  polarimetric  thermal  emission  from  randomly  rough  ocean  surfaces 
using  a  Monte  Carlo  method  has  been  carried  out.  This  study  indicates  that  polarimetric 
information  can  be  useful  in  the  remote  sensing  of  anisotropic  rough  ocean  surfaces.  One 
proposed  application  is  the  remote  sensing  of  wind  direction  over  the  oce  an,  as  the  rough 
surface  created  on  the  ocean  has  an  anisotropic  structure.  Since  the  Ub  parameter  is 
shown  to  be  relatively  insensitive  to  the  polar  angle,  the  permittivity  of  the  medium,  and 
the  cutoff  wavenumber  of  the  surface  spectral  density  function,  a  wind  direction  sensor 
could  feasibly  operate  under  a  wide  variety  of  these  conditions.  Further  research  into  the 
area  and  into  the  use  of  a  more  realistic  ocean  surface  model  will  continue. 
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Figure  1:  Geometry  of  the  problem 


Figure  2:  Convergence  of  polarimetric  brightness  temperature  with  period  of  surfaces  (a) 
UB  (b)  TBh  (c)  TBv  (d)  VB 


Figure  3:  Comparison  of  average  brightness  temperatures  for  20A  and  25A  period  surfaces 
(a)  UB  (b)  Tbh  (c)  TBu  (d)  VB 


Figure  4:  Effect  of  rms  height  on  polarimetric  brightness  temperature  (a)  UB  (b)  TBh  (c) 

tBv  (d)  VB 


Figure  5:  Effect  of  polar  angle  on  polarimetric  brightness  temperature  (a)  UB  (b)  TBh  (c) 
To.  (d)  VB 


Figure  6:  Effect  of  permittivity  on  polarimetric  brightness  temperature  (a)  UB  (b)  TBh  (c) 
TBv  (d)  Vjs 


Figure  7:  Effect  of  power  law  slope  on  polarimetric  brightness  temperature  (a)  UB  (b)  TBh 
(c)  TBv  (d) 


Figure  8:  Effect  of  cutoff  wavenumber  on  polarimetric  brightness  temperature  (a)  UB  (b) 
Tbh  (c)  TBv  (d)  VB 
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Abstract  An  experiment  in  which  the  third  Stokes  parameter  thermal  emission  from 
a  periodic  water  surface  was  measured  is  documented.  This  parameter  is  shown  to  be 
related  to  the  direction  of  periodicity  of  the  periodic  surface  and  to  approach  brightnesses 
of  up  to  30  K  at  X  band  for  the  surface  used  in  the  experiment.  The  surface  actually 
analyzed  was  a  “two-layer”  periodic  surface;  the  theory  of  thermal  emission  from  such 
a  surface  is  derived  and  the  theoretical  results  are  found  to  be  in  good  agreement  with 
the  experimental  measurements.  These  results  further  the  idea  of  using  the  third  Stokes 
parameter  emission  as  an  indicator  of  wind  direction  over  the  ocean. 

1  Introduction 

Recent  theoretical  works  have  suggested  the  potential  of  passive  polarimetry  in  the  remote 
sensing  of  geophysical  media  [1-3].  These  works  indicate  that  the  third  Stokes  parameter, 
UB,  of  the  thermal  emission  may  become  large  for  azimuthally  asymmetric  fields  of  ob¬ 
servation.  In  [3],  values  of  Ub  as  high  as  40  K  were  measured  from  a  periodic  triangular 
soil  surface  at  10  GHz,  and  the  Ub  parameter  was  found  to  indicate  the  azimuthal  angle 
between  the  observation  angle  and  the  direction  of  surface  periodicity.  It  has  thus  been 
proposed  that  these  properties  could  make  passive  polarimetric  measurements  useful  in  the 
remote  sensing  of  wind  direction  over  the  ocean,  since  a  wind  generated  ocean  surface  has 
an  anisotropic  structure.  In  order  to  demonstrate  the  existence  of  the  Ub  parameter  for 
a  water  surface,  and  thus  to  investigate  the  potential  applicability  of  passive  polarimetry 
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to  ocean  remote  sensing,  observations  of  a  periodic  water  surface  were  performed  at  the 
Cold  Regions  Research  and  Engineering  Laboratory  (CRREL),  Hanover,  NH  in  Septem¬ 
ber  of  1992.  This  paper  documents  those  observations  and  the  theory  developed  for  the 
“two-layer”  periodic  surface  actually  analyzed  in  the  experiment. 

2  Theory  of  Polaximetry 

In  passive  polaximetry,  brightness  temperatures  corresponding  to  all  four  modified  Stokes 
parameters  axe  measured.  The  brightness  temperature  Stokes  vector  is  defined  as 

41  (EhE'h)  • 

T  -Lt-L  4  _  J_  m 

B  C  C  u  TjC  2Re(EvE'h)  1 } 

V J  [2Im(EuE^ 

In  the  above  equation,  E h  and  Ev  are  the  emitted  electric  fields  received  from  the  horizontal 
and  vertical  polarization  channels  of  the  radiometer,  T]  is  the  characteristic  impedance, 
and  C  ~  K/X7  with  K  denoting  Boltzmann’s  constant,  A  the  wavelength.  The  first  two 
parameters  of  the  brightness  temperature  Stokes  vector  correspond  to  the  received  powers 
for  horizontal  and  vertical  polarizations,  respectively.  The  third  and  fourth  parameters 
correspond  to  the  complex  correlation  between  the  electric  fields  received  by  the  horizontal 
and  vertical  channels.  We  label  the  four  parameters  Tbhi  TBv,  Ubi  and  Vb  respectively  in 
this  paper. 

It  is  shown  in  [2]  that  the  third  and  fourth  Stokes  parameters  may  be  related  to  the 
brightness  temperature  in  a  45  degree  linearly  polarized  measurement  (TBp)  and  a  right- 
hand-circularly  polarized  measurement  ( TBr )  as  follows: 


UB 

~  2 TBp  -  TBh  —  TBv 

(2) 

Vb 

~  2TBt  -  TBh  —  TBv 

(3) 

Thus,  to  measure  all  four  parameters  of  the  Stokes  vector,  the  brightness  temperatures  in 
horizontal,  vertical,  45  degree  linear,  and  right-hand-circular  polarizations  are  first  mea¬ 
sured,  and  the  above  equations  are  used  to  obtain  Ub  and  Vp. 
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In  the  passive  remote  sensing  of  rough  surfaces,  the  parameter  that  is  actually  of 
interest  is  the  emissivity,  which  relates  the  brightness  temperature  emitted  by  an  object 
to  its  actual  physical  temperature,  under  the  assumption  that  the  object  is  at  a  constant 
physical  temperature,  Tphy,,  and  that  the  emission  from  the  object  is  the  only  source  of 
brightness: 

Tb*  =  ea(9, 4>)Tphyt  (4) 

In  the  above  equation,  the  subscript  a  refers  to  the  polarization  of  the  brightness  tempera¬ 
ture,  6  to  the  polar  observation  angle,  and  <f>  to  the  azimuthal  observation  angle.  Through 
the  principles  of  energy  conservation  and  reciprocity,  Kirchhoff’s  Law  relates  this  emissivity 
to  the  reflectivity  of  the  surface  [4]: 

ea(0,<f>)  =  1  ~ra(0,4>)  .  (5) 

The  reflectivity  ra(0,  <f>)  for  the  given  incident  polarization  a  is  defined  as  the  fraction 
of  the  power  incident  from  direction  (0,  <£)  that  is  rescattered  and  can  be  evaluated  by 
integrating  the  bistatic  scattering  coefficients  76o(^j  4>\  <t>')  over  all  scattering  angles  in 
the  upper  hemisphere  and  summing  the  results  of  both  orthogonal  scattering  polarizations: 

ra(9,4>)  =  JQ  sin  O'  d<j>  7<«(0,  <M',  <£')  (6) 

b 

In  the  expression  of  the  bistatic  scattering  coefficient,  (0,  <f> )  and  (O',  <f>')  represent  the 
incident  and  the  scattered  directions,  respectively,  and  the  subscripts  a  and  b  represent 
the  polarizations  of  the  incident  and  the  scattered  waves,  respectively. 

Thus,  to  calculate  the  fully  polarimetric  emission  vector,  the  bistatic  scattering  co¬ 
efficient  for  each  of  four  polarizations  is  first  calculated  and  integrated  over  the  upper 
hemisphere  to  obtain  the  reflectivity  for  that  particular  polarization.  Multiplication  of  the 
corresponding  emissivity  by  the  physical  temperature  of  the  object  under  view  yields  the 
brightness  temperature  for  this  polarization.  The  fully  polarimetric  brightness  vector  is 
then  calculated  as  described  previously. 
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3  Two-Layer  Periodic  Surface  Theory 

In  order  to  form  a  periodic  water  surface,  a  thin  corrugated  sheet  of  ^berglass  was  placed 
on  top  of  a  flat  water  surface.  After  air  bubbles  trapped  underneath  were  removed,  a 
“two-layer”  periodic  surface  resulted,  as  shown  in  Figure  1.  In  this  figure,  region  zero  is 
the  free  space  region  above  the  periodic  surface,  with  permittivity  eo  and  permeability 
Hq,  region  one  is  the  fiberglass  layer,  with  permittivity  c:  and  permeability  fi i,  and  region 
two  is  the  water,  with  permittivity  e 2  and  permeability  /1 2-  To  examine  the  effects  of  the 
fiberglass  layer  on  the  measured  brightness  temperatures,  a  theoretical  model  for  4hermal 
emission  from  a  two-layer  periodic  surface  was  derived  and  implemented. 

The  problem  of  calculating  thermal  emission  from  a  “single-layer”  periodic  surface  has 
been  solved  by  several  researchers  using  many  different  techniques.  Two  of  the  methods 
used  in  the  single-layer  case  are  extended  for  the  two-layer  case  in  this  analysis.  The  results 
from  these  two  ‘exact’  methods  are  then  compared  to  determine  which  would  apply  best 
to  the  cases  measured  in  the  experiment.  For  both  of  these  methods,  the  calculations  are 
performed  to  determine  the  total  power  reflected  in  the  reciprocal  active  problem;  this 
total  reflectivity  is  then  used  to  obtain  the  emissivity  as  described  previously. 

The  extended  boundary  condition  (EBC)  method  described  by  Chuang  and  Kong  [5] 
for  calculating  scattering  from  a  single-layer  periodic  surface  is  the  first  method  to  be 
applied.  In  this  method,  Huygens’  principle  is  applied  at  the  surface  so  that  the  scattered 
field  can  be  obtained  once  the  surface  fields  are  known.  The  Huygens’  integral  equation  is 
solved  by  expanding  the  unknown  surface  fields  into  Fourier  series  and  forming  a  truncated 
impedance  matrix  which  is  inverted  to  obtain  the  surface  fields.  This  method  is  known  as 
the  extended  boundary  condition  method  because  the  calculations  enforce  the  requirement 
that  the  Huygens’  integral  which  equals  the  field  value  within  the  region  of  interest  must  ue 
zero  outside  that  region.  Although  the  extended  boundary  condition  method  is  exact,  the 
numerical  formulation  becomes  ill  conditioned  if  the  surface  height-to-period  ratio  becomes 
too  large  [10]. 
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To  overcome  this  problem,  a  second  method,  the  method  of  moments  [9],  can  be  used 
to  perform  the  calculations.  In  this  formulation,  the  same  Huygens’  principle  equation  is 
used,  but  the  requirement  that  the  Huygens’  integral .  the  region  not  of  interest  be  zero 
is  not  enforced.  Instead,  the  integral  equation  is  “tested”  at  a  discrete  number  of  points 
along  each  surface  where  the  equation  is  forced  to  hold.  As  the  number  of  testing  points 
is  increased,  the  calculated  scattered  fields  converge  to  the  exact  results.  The  method  of 
moments  can  be  used  in  cases  where  the  extended  boundary  condition  method  fails.  It 
requires,  however,  more  computation  time  than  the  EBC. 

3.1  General  Formulation 

Consider  a  plane  wave  incident  upon  a  two-layer  periodic  surface  described  by  fi(x)  — 
/j(x  +  P ),  /2(x)  —  /2(x  -f  P)  with  /i  representing  the  upper  surface,  /2  the  lower  surface 
and  P  denoting  the  period  of  the  surface  in  the  x  direction  (see  figure  1).  The  electric  field 
of  the  incident  wave  is  given  by 

Ei  =  iiEoe**  (7) 

where  k{  denotes  the  incident  wave  vector  and  is  equal  to  xk^  +  yk^  —  zkzl  and  et-  is  the 
polarization  of  the  electric  field  vector. 

Since  the  structure  is  uniform  in  the  y  direction,  all  the  field  components  in  regions  0,1, 
and  2  will  be  phase  matched  to  the  same  exp  (ifc^y)  dependence.  Thus,  we  can  replace  all 
the  terms  in  Maxwell’s  equations  with  ik ^  as  is  commonly  done  in  waveguide  theory 
(usually  with  the  z  dependence  by  convention.)  Subsequent  equations  will  thus  have  the 
exp  (ikyiy)  dependence  removed. 

Maxwell’s  equations  can  now  be  simplified  so  that  the  x  and  z  components  of  the 
electric  and  magnetic  fields  can  be  expressed  as  functions  of  the  y  components  of  the  fields 
as  follows: 

Ej.{f)  =  ,-2~-7-t fceV, Ejv{f)  +  w/ij-V,  x  U:y{f))  (8) 

Kj  Kyi 

=  ,  ,  *  -  ue,v.  x  £,»(?)]  (9) 

Kj  Kyi 
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where  j  =  0,1,2  signifies  regions  0,1,  or  2,  respectively,  V,  is  the  gradient  operator  that  is 
transverse  to  the  y  direction 

v-  =  ik+ik  (10) 

and  s  =  (x,z)  indicates  the  transverse  components  of  the  fields  for  region  j.  The  y 
components  of  the  electric  and  magnetic  fields  now  satisfy  the  partial  differential  equations 

(v;  +  k?  -  kS)Aju  =  o  (ii) 

where  Ajv  =  Ejv,  Hju.  Since  V,  operates  only  in  a  two-dimensional  space,  the  y  components 
of  the  electric  and  magnetic  fields  satisfy  a  two-dimensional  wave  equation.  The  Green’s 
function  for  such  an  equation  is  given  by 


where  y  =  0, 1,2, 


P.)  =  -  K\) 

(12) 

—  X  exp  [ifc*(z  —  x1)  +  ikjz\z  —  z'|] 

(13) 

kjl  —  yj  kj  kyi 

(14) 

>=i,  =  At  -  ki 

(15) 

p,  =  XX  +  zz 

(16) 

The  above  square  root  functions  are  defined  with  the  branch  cut  on  the  negative  real  axis. 

Appbdng  the  scalar  Green’s  function  form  of  Huygens’  principle  to  the  top  periodic 
surface  and  region  0  gives 

EM.)  -  r  d<T' {G0{p.,-p.)n,  ■  KEo*{P.)  -  Eov{p.)nx  •  V'G0{p.,p.)) 

J  —CO 

__  f  ■E'Ov(Pj)  z  ->  fi(x)  ii?) 

\o  z<h(x)  [  } 

where  the  y  dependence  has  been  suppressed  and 

dcr'fij  —  z  —  x  dx'  (18) 

dx' 
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The  above  surface  integral  is  over  an  infinite  domain.  However,  the  periodic  properties 
of  the  surface  fields  can  be  used  to  reduce  the  integration  domain  to  a  single  period.  The 


surface  fields  have  the  property 


+  xnP )  =  u(p,)e' 


where  u>(pt)  can  be  Ejy(p,),  Hjy(p,)  or  their  normal  derivatives,  n  is  an  integer,  and  P  is 
the  period  of  the  surface. 

Following  Chuang  and  Kong  [5],  a  periodic  Green’s  function  can  be  derived  using  this 
property.  This  Green’s  function  is 

Gjp(p^Ps)  =  2ik~  P  ^  /T~  CXP  ~  **)  +  iki‘Pin  lz  ~  2'IJ  (20) 


where 


kxi  +  Ti^FT 


XI  I  u  P 

ain  =  7 

ki‘ 


"  l  i(!  ~  4 


)1/2  <  <  1 


"  \  i(l  ~  <)l/2  <  >  1  V  ' 

Now  using  the  periodic  Green’s  function  in  the  scalar  Huygens’  principle  formulation 
for  each  of  the  three  regions  yields 
Region  zero 

EM  -  f  da1  {G0P{pM^  •  V.EM)  -  £ov(p>i  •  V.Gop(p.,P.)} 

JP(S  1) 

_  /  £<*(?.)  *  >  /l(*)  m>| 

"10  *  <  Mz)  {23) 


HM  -  [  do' {G0p{p.,p.)h,  •  V.HM)  -  •  V:Gop(p.,pJ> 

JP{Sl) 

I  Eov(ps)  Z  >  fi(z) 

1  o  z  <  /i(x) 
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Region  one 

f  da'  {G,p(p.,K)ni  ■  VMf.)  -  EM)A,  ■  V',G, ,(?„£)}  - 

JP(S  1) 

/  ia' {GlP(p„fi.)h,  ■  v;b„K)  -  s,MVh  ■ 

JP(S2) 


[  0  2  >  /l(z) 

=  ^  El y{p.)  f2(x )  Z  <  fx(x) 

0  z  <  f2(x) 


J  da' {Girtji.,?.)* ,  •  v;fr,„(£)  -  *„(£)*,  •  v:gip(p..p.)}  - 

[  da'  {GMp., £)*«  •  ?',#,„(£)  -  tf„(£)n,  •  V'.GirGvp'.)} 

^P(52) 


=  < 


f  0 
0 


2  >  fl(x) 

/2(a)  <  2  <  /x(r) 

2  <  /2(z) 


Region  two 


Jp(s2)da'{G1PiP.,K)^-  V'.EM)  ~  BMVh  ■  V:Gip(p.,Z)} 


(  E2v(pt) 


2  <  /2O2) 
2  >  fi{*) 


(25) 


(26) 


(27) 


Jp{s2)da'  {G2P{ptJt)n2  ■  V,EM)  -  H2y(p,)h2  ■  V.G2p(pt,?.)} 


_  /  ^2y(P«) 

1° 


2  <  /2(x) 

2  >  /2(z) 


(28) 


where  (5/)  indicates  that  the  integration  is  performed  over  surface  l. 

Once  the  surface  fields  axe  determined,  the  fields  in  regions  0,1, and  2  can  be  obtained 
using  these  relations.  At  present,  there  are  sixteen  unknown  vectors  to  be  obtained:  the  y 
components  of  the  electric  and  magnetic  fields  at  surface  one  for  regions  zero  and  one,  the 
y  components  of  the  electric  and  magnetic  fields  at  surface  two  for  regions  one  and  two, 
and  the  normal  derivatives  of  the  y  component  of  the  electric  and  magnetic  fields  for  each 
of  these  surface-region  combinations.  Electric  and  magnetic  field  boundary  conditions  can 
be  used  at  both  surfaces  to  reduce  the  number  of  unknown  vectors  to  eight. 
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3.2  Boundary  Conditions 

At  both  surfaces,  the  y  component  of  the  electric  and  magnetic  fields  are  tangential,  and 
therefore  must  be  continuous: 


Eoy  -  Ely  z  -  /i(z) 

(29) 

tq 

II 

(30) 

Ely  —  E7y  Z  —  /2(l) 

(31) 

ii 

£ 

(32) 

The  components  of  the  fields  along  the  surface  profile  are 

continuous.  These  relations  can  be  written  as: 

also  tangential  and  must  be 

fii  x  Eo.  =  Hi  X  Ex,  z  =  fx(x) 

(33) 

Hi  x  Ho.  —  hi  x  Hi, 

(34) 

n2  x  Ei.  =  n2  x  E3,  z  =  f2(x) 

(35) 

h2  x  Hi,  =  n2  x  H2, 

(36) 

Note  that  both  the  normal  and  the  s  component  fields  lie  in  the  x  —  z  plane,  so  that 
iij  X  A,  is  a  j  directed  vector  whose  amplitude  corresponds  to  the  tangential  component 
of  the  field  along  the  surface  profile.  Following  Chuang  and  Kong  [5],  the  second  set  of 


boundary  conditions  can  be  written  as: 

y(ni  ■  V.£0y)  =  CoU!  x  V.Hlv  +  c2y(hi  •  V,£ly)  z  —  fi(x)  (37) 

y(hi  ■  V ,Hoy)  =  — dofi-i  x  V,Fik.  +  d2y(hi  •  V ,H\y)  (38) 

y(«2  *  V,Eiv)  =  cifi2  x  V,H2y  +  c3y(n2  •  V  ,E2y)  z  —  f2(x)  (39) 

y(n2  *  V«ifi y)  =  —d\hi  x  "^7  ,E2y  +  dzy[n2  •  V,-fir2lf)  (40) 
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where 


r  k2  k  ■ 

_  l  Ky% 

k\,  u>e0 

_  fiMi 
2  €0  *?. 
kp,  _  ^  kyj 
k\t  u >Ho 

A  -  H iMl 

2  Mo*?. 


ff.2  t  . 

.  „  211*.  _  i  v* 

1  L*I.  °>ei 

.  =  iiMi 

'3  Ci  fc|, 

j  _  *1.  _  -1  *y< 

Ul  j  «/n 

r  _  M**!. 

3  Ml  kj. 


These  eight  boundary  conditions  reduce  the  sixteen  unknown  vectors  to  eight.  The 
remaining  eight  unknowns  are  then  solved  by  expanding  the  surface  fields  in  the  Huygens’ 
integral  equations  into  a  sum  of  unknown  amplitude  basis  functions.  The  functions  used 
in  this  expansion  and  the  equations  solved  differ  for  the  EBC  and  Moment  methods;  each 
method  will  be  treated  separately  beginning  with  the  EBC. . 

3.3  Extended  Boundary  Condition  Method 

H  we  substitute  the  expression  for  the  periodic  Green’s  function  into  the  Huygens’  Principle 
integral  equations  and  factor  out  the  unprimed  exponential  dependence,  the  equations 
reduce  to: 

Region  zero 


0  —  Eyi  —  Ooin 

n 

n 

n  -  u  .  _  V"  J*) 


eiSo»^* 

Vfion 

^  ^  /l(m«) 

(45) 

'T 

> 

E 

V 

n 

(46) 

V,@0n 

*  ^  fl(max) 

(47) 

p'l 

vT^ou 

(48) 

Region  one 


0  ^  ^  ^lln  ryt —  "b  ^  ^  bi2n  r* —  ^  ^  / l(mai) 

n  V  Pin  „  \  Pin 
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EM  —  Z^bun  /7J—  +  X  dim  /7J— 
n  VPln  n  VPln 

a 

f2(max)  Z  fl(min) 

(50) 

1 

0  —  2Lrailn  /73 Z-J  a12r»  /-7\ —  2  /2{m«n) 

n  V  Pin  n  V  Pin 

(51) 

o  =  e  «’  ^  +  £  «BL  ^ 

n  VPln  n  VPln 

...  pttjt  ?.  ...  ^fnP* 

n  VPln  n  VPln 

(52) 

1 

/2(ma«)  ^  ^  ^  /l(min) 

(53) 

» 

n  V'  _(M  c*^,“  ,  V''  _(M  e*ln  „  r 

0  —  /  ,  alln  rn —  +  Z^i  a12n  m —  2  ^  /2(mtn) 

n  VPln  n  VPln 

(54) 

Region  two 

» 

x  e*^2n*^ 

E2yi.Pt}  ^22n  rg —  2  ^  /2(m»n) 

n  VP2n 

(55) 

/  „  ^22n  /7j —  2  ^  ;2(mai) 

n  VPOn 

(56) 

> 

#2„(p.)  =  X)  *22!  /7T-  2  <  /2(mm) 

n  VPln 

(57) 

n  _  V*  -W  €*^*n  ^  x 

”  /  .  a22n  /7j —  2  ^  J2(max) 

n  VP2n 

(58) 

a 

where  the  \z  —  z'\  expression  in  the  periodic  Green’s  function  has  been  replaced  by 

z  —  z' 

when  the  observation  point  is  above  the  surface,  z'  —  z  when  the  observation  point  is 

the  surface.  In  the  above  equations, 

below 

a 

kfn  —  xkj,ajn  ±  zkj,f3jn 

(59) 

> 

”  2ikj,P  Jp(st)  ^  |  "•  ' 

e ~a*m  *•(*')  } 

40-  } 

(60) 

> 

11 

1 - - - .  .....  ...  . . . . . .  . . . . 

1 

. . 

— —  f  d<j' 
2ikj,P  Jp(si) 


■Jfa 


*.  '  V.EjyiK) 


-BM)n, .  v; 


fK. 


and  b^l  and  af£l  have  the  same  form  as  their  electric  field  counterparts  with  Ejy  replaced 
by  Hjy.  We  now  recognize  that  the  scattered  field  in  region  0  is  composed  of  a  discrete  set 
of  plane  waves  traveling  in  the  directions  indicated  by  k These  plane  waves  are  known 
as  the  Floquet  modes;  the  amplitude  of  the  nth  reflected  mode  is  determined  by  60i«- 

For  the  extended  boundary  condition  method,  we  use  the  equations  which  involve  the 
region  of  no  interest  for  each  region.  These  are  equations  (46),(48),(49),(51),(52),(54),(56), 
and  (58).  These  eight  equations  will  uniquely  determine  the  eight  unknown  surface  field 
vectors.  Once  the  surface  fields  are  determined,  the  remaining  equations  are  used  to 
calculate  the  scattered  field  in  each  of  the  three  regions. 

Following  Chuang  and  Kong  [5],  a  Fourier  series  expansion  is  used  for  each  of  the 
unknown  field  quantities  as  follows: 


Surface  one 


Eiy  [/^(z)]  =  53  2o£  ik*x  +  i~E~x 

n  ■* 

27rn 

dcrfi!  •  \/tEiv  [p,(x)j  =  ikudx  Y  2(3* exp  ik^x  +  i——x 

Ely  IMX)1  =  53  2?n  exP  \ikxix  + 

dcrhx  •  V.Hiy  fp,(x)]  =  ikudx  53  26*  exp  ikxix  + 

rt  ■* 


Surface  two 


2nn 

E2V  [/?#(*)]  =  53  2af  exp  ikxix  +  i——x 

n  ■* 

[27TT2- 

ikzix  +  i—p-x 

Eiy  [p,(*)]  =  53  27n  eXP  +  *^T*] 


12 


dcrhi  ■  V,n2y  [^(®)]  =  ik2tdx  ^  25®  exp 

n 


,27rn 

ikxix  4-  i—p~x 


(69) 


Substituting  the  above  expressions  into  the  eight  equations  and  using  the  property 

x  =  S*  (^)  (70) 

leads  to  the  coupled  matrix  equation  of  Appendix  A.  Truncating  the  Fourier  series  of  the 
surface  field  unknowns  results  in  a  finite  matrix  equation  which  can  ue  solved  to  obtain 
the  surface  field  unknowns.  The  upward  going  field  amplitudes  are  then  calculated  to  be 


&oi  =  -CcC^7T  ~  C27T'BifiT  ~  AiaT 

KQ, 

iff  =  doU?aT  -  d2^-%6T  -  XlT 

Kq, 

and  the  downward  going  field  amplitudes  in  region  2  are 

522  =  E2  fiB  +  D2  ocB 

ag>  =  T,iB  + 


(71) 

(72) 

(73) 

(74) 


The  matrices  in  the  above  equations  are  as  described  in  Appendix  A.  The  above  coefficients 
are  the  amplitudes  of  the  y  components  of  the  scattered  and  transmitted  electric  and 
magnetic  fields.  The  £  and  z  components  can  be  calculated  from  equations  (8)  and  (9). 
The  scattered  plane  waves  in  region  zero  propagate  along  the  directions  determined  by 
k^,  while  the  transmitted  plane  waves  in  region  two  travel  along  the  directions  of  k2n. 

Once  the  amplitudes  of  the  scattered  modes  are  calculated,  the  emissivity  of  the  surface 
can  be  obtained  by  integrating  the  total  power  reflected  over  the  upper  hemisphere  and 
then  applying  Kirchhoff’s  law  to  obtain  the  emissivity  corresponding  to  the  polarization 
of  the  incident  field.  For  the  nth  propagating  scattered  Floquet  mode,  the  reflectivity  is 
given  by  power  transmitted  in  the  z  direction  divided  by  the  incident  power  and  is  found 
in  [5]  to  be: 


Ifeoinl2  +  hofeffj2 

|ooro|2  +  |7?oa-oioi2 
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where  the  aoio  and  ct^io  ^  35  defined  in  Appendix  A.  The  brightness  temperature  of  the 
surface  is  then  given  by 

TBa  =  W  l-£rn)  (76) 

n 

where  the  sum  is  over  the  propagating  modes  reflected  for  an  incident  wave  polarized  in 
the  a  direction. 


3.4  Method  of  Moments 


A  second  method  for  solving  the  Huygens’  integral  equations  for  a  two-layer  periodic  surface 
is  the  method  of  moments.  The  use  of  the  method  of  moments  for  computing  scattering 
from  a  single-layer  periodic  surface  has  been  presented  by  Veysoglu  et  al  [2];  these  methods 
will  be  extended  for  the  two-layer  periodic  surface  in  this  section. 

We  begin  with  the  integral  equations  presented  in  the  general  formulation  section  and 
the  periodic  Green’s  function.  We  will  also  use  a  second  form  of  the  periodic  Green’s 
function  in  the  calculations 

Gjp(p,p)  =  ~  £  Eil\kj,y/[z  -  (r'+mp)]2  +  (z  -  z')2)exp (ik0xmp)  (77) 

m~— *oo 

in  addition  to  the  sum  of  exponentials  form  used  in  the  EBC  method.  The  eight  surface 
field  unknowns  in  the  integral  equations  are  now  expanded  into  sums  of  unknown  amplitude 
“pulse”  functions  as  follows: 

Surface  one 


Surface  two 


Eiv(r) 

=  £o nPn(f) 

n 

(78) 

fii  *  V,Eiy(r) 

= 

n 

(79) 

Bi  v(r) 

=  ]£cnPn(f) 

n 

(80) 

fix  •  V,fflv(r) 

=  EW) 

n 

(81) 

E2y(f) 

=  X>n-Pn(f) 

(82) 

« 
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(83) 

(84) 

(85) 


n2  ■  V.E„(f)  =  2/»P-(f) 

n 

H2y{r)  =  Yj 9nPn(f ) 

n 

n2‘VtH2v(r)  =  £>„Pn(r) 


where 


W) 


1  £?-.*)£  <  x  <  nE. 

1  N  ^  x  ^  N 


=  {1 

\  0  c 


otherwise 


and  jV  is  the  total  number  of  basis  (pulse)  functions  used. 

Instead  of  solving  the  equations  with  right-hand-sides  of  zero,  as  was  done  in  the  EBC 
method,  the  method  of  moments  uses  the  equations  with  right-hand-sides  corresponding 
to  the  fields  in  the  region  of  interest.  These  six  equations  are  “tested”  at  a  discrete  set  of 
points  along  the  surface  upon  which  the  integral  is  performed  where  the  integral  equations 
are  forced  to  hold.  The  missing  two  equations  needed  to  match  the  eight  unknown  functions 
are  obtained  from  testing  the  integral  equations  for  region  one  on  both  surfaces  one  and 
two.  The  testing  points  are  chosen  to  lie  in  the  center  of  the  pulse  basis  functions  described 
above.  The  integral  equations  now  reduce  to  the  matrix  equation  given  in  the  Appendix 


Note  that  this  method  will  require  the  integration  of  the  singularity  in  the  periodic 
Green’s  function  when  the  testing  point  and  integration  range  overlap.  For  these  “self” 
terms  of  the  matrices,  an  asymptotic  expression  for  the  Hankel  function  of  a  small  argument 
is  used  and  analytically  integrated.  The  integral  expression  proposed  by  Veysoglu  et  al  [2] 
is  also  used  to  speed  up  the  convergence  of  the  evaluation  of  the  Hankel  function  sum  in 
the  periodic  Green’s  function.  After  solving  the  matrix  equation  to  determine  the  surface 
fields,  the  coefficients  of  the  reflected  and  transmitted  Floquet  modes  can  be  calculated. 
If  the  reflected  fields  in  region  zero  are  written  as 

K.  =  E‘««*v' 
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(87) 

(88) 


and  the  transmitted  fields  in  region  two  as 


K 


n 

H\v  = 


the  coeffidents  are  given  in  terms  of  the  surface  fields  by 


6n  = 


2ik0,P  JP(Si) 


[  dS 

JP(Sl) 


Z-'Kn'P 


-ni  ■  VaE0y(p) 


-Eoyififl!  •  V,- 


bW  =  — l —  f  £s 
n  2ik0,P  Jp(si) 


e  *^0n 

~x 


fil  •  V,-GToy(/>) 


-Hoy(p)h I  •  V, 


~~K7 


=  _!_/ 
2ik2,P  Jf 


dS 


P(S2) 


•  V ,E2v(p) 


-E2v(p)fi2  •  V, 


e-**a*'^ 


= 


2ikz,P  JP(S2) 


L 


dS 


n2  *  V.ff2v(p) 


~32v(p)h2  •  V, 


(89) 

(90) 


(91) 


(92) 


(93) 


(94) 


where  is  defined  as  in  the  EBC  section. 

As  with  the  EBC,  once  the  mode  amplitudes  are  determined,  the  emissivity  can  be 
found  by  integrating  the  total  power  reflected  over  the  upper  hemisphere  and  then  applying 
Kirchhoff ’s  Law.  Note  that  the  method  of  moments  does  not  suffer  from  the  problems  of  the 
EBC  when  the  height-to-period  ratio  becomes  too  large.  A  larger  number  of  basis  functions 
is  required  and  the  calculations  take  longer  to  perform,  but  the  method  ultimately  yields 
an  accurate  result. 
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3.5  Comparison  of  Methods 

Both  the  EBC  and  MOM  methods  were  used  to  perform  a  theoretical  analysis  of  the  surface 
in  the  experiment.  Figure  2  is  a  plot  of  the  four  polarimetric  brightness  temperatures 
calculated  by  both  the  EBC  and  MOM  at  polar  angle  20°.  Also  shown  in  these  plots 
are  the  corresponding  brightness  temperatures  for  a  periodic  water  surface  without  the 
fiberglass  layer  on  top.  These  calculations  axe  for  the  experimental  parameters  described 
in  the  next  section,  namely  a  6.8  cm  period  quasi-sinusoidal  surface  (described  in  the  next 
section)  with  a  peak-to-peak  height  of  1.4  cm  and  thickness  1.75  mm,  dielectric  constant 
3.1  +  z0.05  of  the  fiberglass  layer,  and  dielectric  constant  64  + 129  for  the  water,  assuming 
a  physical  temperature  of  300  K.  The  EBC  and  MOM  results  are  seen  to  agree,  indicating 
that  the  surface  slope  limitations  of  the  EBC  are  not  being  exceeded  and  that  the  numerical 
resolution  of  the  method  of  moments  code  is  adequate.  The  MOM  code,  however,  was  found 
to  require  approximately  ten  times  more  CPU  time  than  the  EBC,  primarily  due  to  the 
summation  of  the  slowly  convergent  periodic  Green’s  function.  The  more  efficient  EBC 
method  was  thus  used  for  the  rest  of  the  theoretical  calculations  presented  in  this  paper. 
It  is  also  seen  that  the  effect  of  the  layer  is  to  increase  Tbh ,  Tbv,  and  Tbp  by  approximate1'/ 
30 K  over  the  single-layer  brightnesses,  but  the  effect  on  Ub  is  small  due  to  cancellation 
of  this  increase  when  the  brightnesses  in  the  three  polarizations  are  used  to  compute  Ub • 
Finally,  the  characteristic  “jumps”  in  brightness  temperature  for  an  infinite  periodic  surface 
due  to  the  transition  of  a  Floquet  mode  from  propagating  to  non-propagating  are  seen  in 
the  Tbh  curve.  Note  that  the  exact  shape  of  the  curve  may  not  be  the  same  as  that  of 
Figure  2,  as  the  curve  was  linearly  interpolated  between  the  calculated  points  every  15 
degrees  in  azimuth. 

4  Experiment  Setup 

To  form  the  periodic  water  surface,  a  corrugated  sheet  of  fiberglass  1.32  m  x  1.22  m  x 
1.75  mm  thick  was  placed  on  top  of  a  1.83  m  side  square  pool  of  fresh  water  of  0.23  m 
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depth.  The  corrugation  of  the  sheet  had  a  sinusoidal  profile  on  its  lower  surface,  and  had 
a  uniform  thickness  normal  to  this  surface,  which  defined  the  upper  profile.  This  upper 
surface  is  thus  a  slight  deviation  from  a  true  sinusoid;  we  will  describe  this  surface  as  quasi- 
sinusoidal.  The  height  and  period  of  the  sinusoidal  lower  profile  were  1.4  cm  peak-to-peak 
and  6.8  cm  respectively.  To  prevent  water  from  flowing  over  the  edges,  a  styrofoam  border 
approximately  3  cm  high  was  glued  around  the  surface. 

The  temperature  of  the  water  was  measured  around  the  surface  in  the  experiment  and 
found  to  be  22.5°  ±  1.5°  C  throughout  the  experiment.  As  fresh  water  was  used,  this 
corresponds  to  an  average  dielectric  constant  of  64  +  i29  at  10  GHz  over  this  temperature 
range  [7].  The  dielectric  constant  of  the  fiberglass  layer  was  measured  using  a  network 
analyzer  technique  [6]  from  12  to  18  GHz  and  found  to  be  3.1  ±  0.1  +  i(0.05  ±  0.05)  over 
this  entire  band.  An  accurate  calibration  kit  was  not  available  for  the  network  analyzer 
at  10  GHz;  however,  due  to  the  constant  nature  of  the  dielectric  constant  Horn  12  to  18 
GHz,  it  was  assumed  that  the  dielectric  constant  was  the  same  at  10  GHz. 

The  radiometer  used  in  the  experiment  is  the  same  as  that  used  in  [3].  It  operates 
at  10  GHz  with  a  200-MHz  RF  bandwidth,  —3  mV/K  sensitivity,  ±1  K  precision,  and 
a  pyramidal  (standard  gain)  horn  of  30°  beamwidth.  The  radiometer  was  mounted  on  a 
tripod  at  an  elevation  of  1.7  m  above  the  pool  surface  and  directed  toward  the  surface 
along  a  direction  determined  by  the  azimuthal  angle,  <f>,  and  the  polar  angle,  0  (see  Figure 
1).  The  radiometer  was  calibrated  by  viewing  an  absorber  whose  physical  temperature 
was  known  for  the  hot  load,  and  by  viewing  the  sky  at  a  set  of  angles  for  the  cold  load. 
Since  the  measured  sky  voltage  can  be  expressed  as  a  function  of  polar  angle,  the  points 
measured  can  be  extrapolated  to  obtain  a  radiometer  voltage  corresponding  to  brightness 
temperature  3 K.  This  voltage  is  then  used  as  the  cold  reference  for  the  calibration. 

Observations  of  the  periodic  surface  were  made  as  a  function  of  azimuthal  angle,  which 
was  varied  by  rotating  the  periodic  surface  in  the  pool  so  that  the  background  noise  contri¬ 
bution  from  the  antenna  sidelobes  would  remain  constant  through  a  set  of  measurements. 
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The  horn  of  this  single  polarization  radiometer  was  rotated  to  prodm  measurements  of 
Tbh,  2bu,  and  Tbp •  Measurements  of  Tbt  were  not  possible  with  this  radiometer;  however, 
the  Vb  brightness  was  not  predicted  to  be  significant. 

5  Estimation  of  Noise  Contribution 

The  contribution  of  background  noise  to  the  measured  antenna  temperature  for  the  periodic 
surface  can  be  large:  noise  from  the  antenna  sidelobes  and  from  reflections  of  sky  bright¬ 
ness  off  the  surface  both  contribute  to  error  in  the  measurements.  Since  this  background 
noise  could  not  be  easily  eliminated  from  the  experiment,  estimates  of  tbe  background 
contribution  were  made  so  that  the  noise  could  later  be  removed.  The  contribution  of  the 
background  noise  was  estimated  by  measuring  the  brightness  temperature  of  a  flat  water 
surface  alone  and  by  measuring  the  brightness  temperature  of  reflector  in  the  pool  which 
was  the  same  size  as  the  periodic  surface.  The  theoretical  results  for  both  of  these  cases 
are  known;  thus,  these  measurements  form  a  sort  of  two-point  calibration  for  the  periodic 
surface  when  it  is  in  the  pool. 

Consider  the  case  of  the  observation  of  the  water  in  the  pool  alone.  The  measured 
brightness  temperature,  consists  of  contributions  from  the  water  sur'ice,  Tjiat, 

from  the  ground  beside  the  pool  and  the  sides  of  the  pool,  T,itu,  and  from  reflections  of 
the  sky  temperature,  T,ky,  off  the  pool.  The  contribution  of  the  reflections  of  the  sky 
temperature  from  the  ground  beside  the  pool  is  neglected.  If  we  assume  that  there  is  no 
variation  in  these  temperatures  with  the  incidence  angle,  then  the  relationship  between 
these  variables  can  be  written  as: 

TfiaiM  =  —  /,)  +  T,kv(  1  -  /,)( 1  -  Tflat/TPwat )  -j-  T,idcf,  (95) 

where  f,  represents  the  fraction  of  the  antenna  pattern  that  does  not  view  the  pool  and 
Tp^at  is  the  physical  temperature  of  the  water  in  the  pool.  Measured  values  of  f,  ranged 
from  3  to  18  percent  through  the  experiment.  Both  T,kv  and  T,i,u  were  obtained  from  the 


experiment  at  angles  corresponding  to  the  specular  reflection  for  the  sky  temperature  and 
direct  incidence  for  the  ground  temperature. 

Next  consider  the  case  of  a  reflector  in  the  pool  the  same  size  and  at  the  same  loca¬ 
tion  as  the  periodic  surface.  The  brightness  temperature  measured,  Tre/tM,  consists  of  the 
reflection  of  the  sky  temperature  over  the  reflector  surface,  the  flat  water  surface  temper¬ 
ature  over  the  remainder  of  the  pool  and  the  reflection  of  the  sky  temperature  from  this 
portion  of  the  surface,  and  the  ground  contribution.  This  can  be  expressed  as: 

TrcflM  —  T.kyfr  +  T,ky(l  ~  fr  ~  f »)( 1  ~  Tflat/Tpwat ) 

—  /.  —  /«)  +  T,i<uf,  (96) 

=  TfutM  +  Tjiaifr(T,ky/Tpwat-  1)  (97) 

where  fr  is  the  fraction  of  the  antenna  pattern  that  views  the  reflector  surface,  which  was 
found  to  range  from  50  to  74  percent  through  the  experiment.  A  further  assumption  that 
the  physical  temperature  of  the  water  and  the  background  sky  temperature  are  the  same 
in  the  flat  surface  and  reflector  measurements  is  made  in  the  second  equality. 

Finally,  consider  the  case  of  the  periodic  surface  in  the  pool.  The  measured  tempera¬ 
ture,  T.ut}M,  consists  of  contribution  from  the  periodic  surface,  the  flat  water  in  the  pool, 
sky  reflections  off  the  pool  area,  and  the  side  contribution.  The  reflections  of  sky  tem¬ 
perature  off  the  periodic  surface  are  determined  by  its  Floquet  modes;  however,  since  the 
sky  temperature  is  assumed  uniform,  this  reduces  to  the  same  case  as  that  of  a  specular 
reflector.  The  measured  temperature  can  be  expressed  as: 

Tiur/A/  —  T,urffr  4-  T,kyfr(l  ~  T,urf/Tpwat )  +  T flat{  1  ~  ft  ~  ft) 

+T.kv(  1  -  Tflat/TPwat)(  1  -  fT  -  f.)  +  T.ukf.  (98) 

=  'TrefiM  +  Trurffr(  1  —  T,^y/Tp  wat  )  (99) 


These  three  equations  can  be  used  to  solve  for  T,urf ,  the  periodic  surface  brightness 


temperature  averaged  over  the  fraction  of  the  antenna  pattern  on  the  surface: 

T,urj  —  ~ - (T,urfM  -  Tre}m)  (100) 

This  brightness  temperature  has  the  effects  of  both  the  reflections  of  the  sky  temperature 
off  the  surface  and  the  antenna  field  of  view  outside  of  the  periodic  surface  removed. 

This  calibration  procedure  increases  the  effect  of  the  radiometer  measurement  uncer¬ 
tainty  in  the  final  calculated  Tturj  to  approximately  3  K.  The  assumptions  made  in  the 
above  equations  also  contribute  to  the  error  in  the  experiment.  A  theoretical  study  in 
which  the  brightness  contributions  of  the  water  around  the  surface  and  the  sky  tempera¬ 
ture  were  estimated  by  integrating  the  angularly  varying  T/jat  and  T,kv  over  a  theoretical 
antenna  pattern  indicated  that  the  use  of  the  specular  Tfiat  and  Ttkv  values  were  valid  to 
within  2  K.  A  bias  in  both  of  these  values  would  produce  a  stematic  error  of  less  than 

5  K.  One  remaining  possible  source  of  error  is  the  possibility  of  the  close-to-the-horizon 
Floquet  modes  of  the  periodic  surface  imaging  something  with  a  higher  brightness  temper¬ 
ature  them  the  sky.  Note  that  the  Ub  parameter  is  fairly  insensitive  to  any  systematic  error 
in  the  calibration  used.  Since  systematic  error  will  tend  to  affect  all  three  polarizations  in 
a  similar  manner,  the  effects  tend  to  cancel  out,  leaving  Ub  unaffected. 

The  results  presented  in  Section  6  for  the  measurements  of  periodic  surface  brightness 
temperatures  were  obtained  using  the  above  technique. 

6  Experimental  Results 

Figure  3  is  a  plot  of  the  measured  brightness  temperatures  before  and  after  the  noise  re¬ 
moval  calibration  for  8  —  20°,  normalized  to  a  constant  physical  temperature  of  300  K. 
Plotted  are  the  measured  brightness  temperatures  Tbh,Tbv ,  Ub,  and  Tbp.  The  Ub  bright¬ 
ness  temperature  is  obtained  from  Ub  —  Tsh  +  Tbv  —  2Tbp-  This  formula  is  opposite  in  sign 
to  that  previously  given  due  to  the  fact  that  the  radiometer  polarization  basis  was  actually 
rotated  to  —45°  in  the  experiment.  It  is  observed  from  the  experiment  that,  in  general, 
the  Ub  value  approaches  a  maximum  absolute  value  at  <$>  —  —45°  and  approaches  zero  at 
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(ft  =  —90°  and  <f>  =  0°.  This  shows  that  the  Ub  parameter  is  sensitive  to  the  azimuthal 
direction  of  a  periodic  water  surface.  It  is  also  observed  that  the  noise  removal  calibration 
has  its  largest  effects  near  <fi  =  —45°,  at  which  the  lowest  amount  of  the  radiometer  beam 
pattern  was  on  the  surface. 

Figure  4  plots  the  measured  brightness  temperatures  at  $  =  30°  before  and  after  the 
noise  removal  calibration.  The  same  trends  in  the  Ub  temperature  axe  observed  as  in 
Figure  3.  Note  that  the  noise  removal  technique  has  a  larger  effect  on  the  data  for  this 
case,  due  to  the  fact  that  less  of  the  antenna  pattern  was  on  the  surface  for  the  higher 
polar  angle. 

7  Comparison  with  Plane  Wave  Calculations 

Figures  5  and  6  are  plots  of  the  measurement  results  and  the  theoretically  predicted  bright¬ 
ness  temperatures  for  polar  angles  20°  and  30°,  respectively,  normalized  to  a  constant  phys¬ 
ical  temperature  of  300  K.  The  theoretical  values  plotted  were  obtained  using  the  modified 
version  of  the  Extended  Boundary  Condition  method.  The  theoretical  and  experimental 
results  are  seen  to  have  the  same  general  trends  in  all  both  of  these  figures.  However,  the 
abrupt  changes  in  the  theoretical  curves  due  to  the  transition  of  one  Floquet  mode  of  the 
periodic  surface  from  propagating  to  non-propagating  are  not  seen  in  the  experimental 
results.  The  primary  reason  for  this  discrepancy  is  considered  in  the  next  section. 

8  Antenna  Pattern  Effects 

One  source  of  the  differences  between  the  theoretical  and  experimental  results  presented 
in  the  last  section  is  the  effect  of  the  radiometer  antenna  pattern.  The  theoretical  results 
presented  in  the  last  section  assumed  a  plane  wave  incidence  or,  for  the  passive  case,  an 
antenna  pattern  which  viewed  only  one  angle  (a  delta  function  in  theta  and  phi  at  the 
observation  angle.)  Since  the  antenna  used  in  the  experiment  actually  samples  a  range  of 
angles  over  its  pattern,  this  effect  should  be  taken  into  account  in  order  to  more  accurately 


22 


model  the  experiment.  To  take  the  antenna  pattern  effects  into  account,  the  observed 
antenna  temperature,  Tx,  can  be  written  as: 

A  1 1  g(o,  f)<in 

where  G(8 ,  4>)  is  the  power  gain  pattern  of  the  antenna  and  the  integrals  are  performed 
over  all  space.  This  expression  assumes  that  the  emission  from  different  directions  is 
uncorrelated,  since  the  brightness  temperatures  (powers)  are  summed  directly. 

The  gain  pattern  used  for  the  calculations  was  obtained  from  [8],  using  the  pattern  for 
a  pyramidal  horn  with  the  measured  antenna  dimensions  of  b  =  5.7  cm  by  a  =  7  cm  for  the 
aperture  dimension  and  Ip  —  10  cm  and  I3  —  8.9  cm  for  the  corresponding  flare  lengths 
of  the  horn. 

Figures  7  and  8  are  plots  of  the  measurement  results  and  the  theoretically  predicted 
values  for  the  experiment  at  polar  angles  20°  and  30°,  respectively.  The  theoretical  values 
plotted  were  obtained  using  the  modified  version  of  the  Extended  Boundary  Condition 
method.  These  theoretical  results  were  also  averaged  noncoherently  over  81  points  in  the 
three-dimensional  calculated  antenna  pattern  for  the  pyramidal  horn  used  in  the  experi¬ 
ment,  which  tends  to  smooth  out  the  abrupt  changes  obtained  for  plane  wave  incidence  in 
Section  7. 

Good  agreement  is  observed  between  the  theory  and  experiment,  indicating  that  the 
background  noise  removal  procedure  was  reasonable.  However,  there  is  a  slight  bias  of 
the  experimental  results  above  the  theoretical  predictions  on  the  average.  One  possible 
source  of  this  error  is  the  dielectric  constant  value  for  the  fiberglass  surface  used  in  the 
experiment,  which  could  be  slightly  higher  at  10  GHz  than  the  measured  value  at  12.4 
GHz.  An  increase  in  the  real  part  of  this  value  to  3.2  indicates  that  this  could  increase 
the  theoretical  predictions  by  approximately  5  K  or  0.017  emissivity.  A  second  possible 
source  of  the  error  is  the  use  of  the  specular  values  for  T//at  and  T,kv  as  discussed  in 
section  5.  This  could  also  contribute  approximately  5  K  or  0.017  emissivity  variation. 
Finally,  the  imaging  of  a  brightness  source  other  than  the  sky  by  the  Floquet  modes  of  the 
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periodic  surface  would  tend  to  increase  the  experimental  brightness  temperatures  above 
the  theoretical  predictions.  However,  the  area  in  which  the  experiment  was  performed  was 
clear  to  about  ten  degrees  above  the  horizon.  Only  the  azimuthal  angles  —90,  —15,  0,  and 
15  are  theoretically  predicted  to  have  modes  this  close  to  the  horizon. 

9  Conclusion 

We  have  demonstrated  in  this  paper  that  the  Ub  parameter  exists  for  a  two-layer  periodic 
water  surface  and  can  approach  brightnesses  exceeding  30  K  at  X  band.  We  have  also 
shown  that  the  U b  parameter  is  fairly  insensitive  to  the  effects  of  the  fiberglass  layer  and  the 
measurement  uncertainties  in  the  experiment  as  an  indicator  of  surface  azimuthal  direction. 
Although  the  periodic  water  surface  is  an  extremely  simplified  model  of  an  actual  ocean 
surface,  this  experiment  further  strengthens  the  idea  of  using  passive  polarimetry  to  infer 
wind  direction  over  the  ocean.  Further  research  into  this  area  and  into  other  applications 
cf  passive  polarimetric  remote  sensing  will  continue. 
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Appendix  A:  EBC  Coupled  Matrix  Equation 


The  matrix  equation  of  the  EBC  can  be  written  as: 


where 


1  r  e-**i«*w  2n7r 
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aAEt), 


(108) 


where  ( Sj )  indicates  that  the  integration  is  over  surface  j, 

UOln  =  ^riO^/^O-EoCfit  ‘  y) 


(109) 


4ln  =  GnO)f(3o~(h  X  Bi  -y) 


(110) 


The  above  5oi  elements  are  determined  by  the  y  component  of  the  incident  electric  and 
magnetic  fields.  Truncating  the  Fourier  series  of  the  surface  field  unknowns  results  in  a 
finite  matrix  which  can  be  inverted  to  obtain  the  surface  field  unknowns. 


The  integrals  of  the  above  matrix  were  evaluated  exactly  over  a  linearly  interpolated 
sampled  profile  of  the  sinusoidal  and  quasi-sinusoidal  surfaces. 
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Appendix  B:  MoM  Matrix  Equation 


The  matrix  equation  for  the  Method  of  Moments  can  be  written  as: 
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with  the  testing  points  on  surface  one,  integration  on  surface  two, 
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with  the  testing  points  on  surface  two,  integration  on  surface  one,  and 
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with  the  testing  points  and  integration  on  surface  two.  The  above  integrals  were  approxi¬ 
mated  by  a  single  point  rectangular  area  at  the  testing  point.  This  approximation  becomes 
more  accurate  as  the  number  of  testing  points  is  increased.  In  the  results  presented,  150 
basis  functions  were  used,  which  corresponds  to  66  basis  functions  per  wavelength  in  the 
free  space  region.  This  large  number  of  basis  functions  was  needed  due  to  the  much  shorter 
wavelength  in  the  high  permittivity  water  region. 


10  Figure  Captions 


Figure  1:  Geometry  of  a  “two-layer”  periodic  surface 


Figure  2:  Comparison  of  EBC  and  MOM  brightness  temperatures  at  polar  angle  20°  versus 
azimuthal  angle  (a)  UB  (b)  TBh  (c)  TBv  (d)  VB 


Figure  3:  Measured  brightness  temperatures  before  and  after  noise  removal  calibration  at 
polar  angle  20°  versus  azimuthal  angle  (a)  UB  (b)  Tbk  (c)  TBv  (d)  VB 


Figure  4:  Measured  brightness  temperatures  before  and  after  noise  removal  calibration  at 
polar  angle  30°  versus  azimuthal  angle  (a)  UB  (b)  TBh  (c)  TBv  (d)  VB 


Figure  5:  Comparison  of  theoretical  and  experimental  brightness  temperatures  at  polar 
angle  20°  versus  azimuthal  angle  (a)  UB  (b)  TBh  (c)  TBv  (d)  VB 


Figure  6:  Comparison  of  theoretical  and  experimental  brightness  temperatures  at  polar 
angle  30°  versus  azimuthal  angle  (a)  UB  (b)  TBh  (c)  TBv  (d)  VB 


Figure  7:  Comparison  of  beam  averaged  theoretical  and  experimental  results  for  polar 
angle  20  degrees  (a)  Ub  (b)  TBh  (c)  TBv  (d)  TBp 


Figure  8:  Comparison  of  beam  averaged  theoretical  and  experimental  iesults  for  polar 
angle  30  degrees  (a)  UB  (b)  TBh  (c)  TBu  (d)  TBp 
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